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Abstract 

We use the duality between color and kinematics to simplify the construction of the complete 
four- loop four-point amplitude of M = 4 super- Yang-Mills theory, including the nonplanar contri- 
butions. The duality completely determines the amplitude's integrand in terms of just two planar 
graphs. The existence of a manifestly dual gauge-theory amplitude trivializes the construction 
of the corresponding N = 8 supergravity integrand, whose graph numerators are double copies 
(squares) of the N = 4 super- Yang-Mills numerators. The success of this procedure provides fur- 
ther nontrivial evidence that the duality and double-copy properties hold at loop level. The new 
form of the four-loop four-point supergravity amplitude makes manifest the same ultraviolet power 
counting as the corresponding N = 4 super- Yang-Mills amplitude. We determine the amplitude's 
ultraviolet pole in the critical dimension of D = 11/2, the same dimension as for J\f = 4 super- Yang- 
Mills theory. Strikingly, exactly the same combination of vacuum integrals (after simplification) 
describes the ultraviolet divergence of J\f = 8 supergravity as the subleading-in-1/./V,? single-trace 
divergence in N = 4 super- Yang-Mills theory. 

PACS numbers: 04.65. +e, 11.15.Bt, 11.30.Pb, 11.55.Bq 
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I. INTRODUCTION 



The past few years have brought remarkable advances in understanding scattering am- 
plitudes in the maximally supersymmetric Af = 4 super- Yang-Mills (sYM) theory [1] in the 
planar limit of a large number of colors. It may soon be possible to completely determine 
all planar scattering amplitudes in this theory, for all values of the coupling, going far be- 
yond the (now thoroughly understood) cases of four and five external gluons [2]. Much of 
this progress has been surveyed recently [3]. Planar scattering amplitudes exhibit a new 
symmetry known as dual conformal symmetry [4, 5], which severely restricts their structure. 
Together with supersymmetry and (position space) conformal symmetry, dual conformal 
invariance gives rise to a Yangian [6] — an algebraic structure common in integrable models. 
Indeed, it is widely believed that several aspects of the planar sector of Af = 4 sYM theory 
are controlled by an integrable model (see e.g. ref. [7]). 

In contrast, much less is known about the nonplanar sector — the subject of the present 
paper. Consider Af = 4 sYM theory for the gauge group SU(N C ). In the limit N c — > oo, the 
nonplanar, or subleading-color, contributions are suppressed by powers of 1/N C . Once one 
takes into account these corrections, for finite N c , the scattering amplitudes no longer appear 
to possess dual conformal symmetry, nor do they demonstrate any obvious integrability 
properties. 

Understanding the subleading-color terms is critical to a complete description of the 
behavior of gauge theories. For example, many types of color correlations are suppressed 
in the large- N c limit. Furthermore, the information provided by the full-color expression 
for Af = 4 sYM amplitudes, expressed in terms of their loop-momentum integrands, can be 
used to construct corresponding amplitudes [8-11] in Af = 8 supergravity [12]. From each 
set of amplitudes one can extract information about ultraviolet divergences in the respective 
theory. 

The ultraviolet (UV) properties of Af = 8 supergravity have been the focus of intense 
investigation. There have been several recent reviews of the situation [13]. Long ago, 
an Af = 8 supersymmetric local counterterm at three loops in D = 4 was proposed [14- 
18]. An explicit computation of the three-loop four-graviton amplitude first revealed that 
the counterterm has a vanishing coefficient [9]. Subsequently it was realized [19] that this 
counterterm is forbidden in D = 4 by the £7(7) duality symmetry [12]. Other analyses have 
extended the finiteness constraints from £7(7) and linearized supersymmetry, such that the 
first potential divergence in D = 4 is now at seven loops [20-23]. Finiteness until this loop 
order happens to agree with an earlier naive power-counting, based on the assumption of an 
off-shell Af = 8 superspace [24]. A potential seven- loop divergence is also suggested by other 
approaches, including an analysis of string theory dualities [25], a first-quantized world-line 
approach [26], and light-cone supergraphs [27]. However, it has also been argued that the 
theory may remain finite beyond seven loops [28]. 

In this paper, we will show how a conjectured duality between color and kinematics [29, 30] 
provides a powerful method for computing subleading-color terms in Af = 4 sYM amplitudes, 
in a way that makes the construction of the corresponding Af = 8 supergravity amplitudes 
extremely simple. Also, the Af = 8 result is expressed in a form that makes manifest the 
true ultraviolet behavior of the amplitude (when continued to higher space-time dimension 
D). Thus this method provides unprecedented access to the precise coefficients of potential 
counterterms in Af = 8 supergravity, as well as in its higher- dimensional versions. It may 
eventually offer a means for settling the question of whether additional UV cancellations exist 
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in M = 8 supergravity, beyond the known or expected ones. Perhaps even more importantly, 
the method gives a means for constructing complete amplitudes, allowing for detailed studies 
of their symmetries and properties. 

A key point is that when the color-kinematics duality holds manifestly, it locks the non- 
planar contributions to the planar ones. The nonplanar contributions are essential for eval- 
uating gravity amplitudes, because in gravity theories no separation exists between planar 
and nonplanar contributions. This duality allows one to efficiently export information from 
the planar sector, e.g. that provided by dual conformal symmetry, to the much more intricate 
nonplanar sector. 

A second key point is the claim [30] that if a duality-respecting representation of M = 4 
sYM amplitudes can be found, then the loop-momentum integrands of the corresponding 
M = 8 supergravity amplitudes can be obtained simply by taking the graphs of M = 4 sYM 
theory, dropping the color factors and squaring their kinematic numerators. This double- 
copy property is a loop-level generalization of the corresponding tree-level property [29], 
equivalent to the Kawai-Lewellen-Tye (KLT) relations between gravity and gauge-theory 
amplitudes [31]. Using the color- kinematics duality, followed by the double-copy property, 
advances in constructing integrands for the planar sector of gauge theory can be carried over 
to the nonplanar sector, and then on to gravity. The color-kinematic duality and the grav- 
ity double-copy property do not appear to require supersymmetry, although amplitudes in 
supersymmetric theories are generally much simpler to work with than non-supersymmetric 
amplitudes. Another important aspect of the duality is that it appears to hold in any 
dimension, thus making it compatible with dimensional regularization. 

In this paper we will exploit the color-kinematic duality to construct the complete four- 
loop four-point amplitudes of M = 4 sYM theory and M = 8 supergravity. Both amplitudes 
were constructed previously by us [10, 11]; however, the present construction is considerably 
more efficient, and makes various properties of the amplitudes more manifest. The color- 
kinematic duality relations allow us to express the four-loop loop-momentum integrands, 
for 83 different cubic (trivalent) graphs, as functionals of the integrands of just two planar 
graphs. For the one-, two- and three-loop four-point [30, 32], and the one- and two-loop 
five-point cases [33], the duality is even more restrictive: a single planar graph suffices 
to determine all the others. As it is becoming increasingly simple to construct planar 
amplitudes, a particularly attractive aspect of using the color-kinematic duality is that 
it determines nonplanar contributions from planar ones. Perhaps even more remarkable, 
in terms of measuring the redundancy found in local gauge-theory scattering amplitudes, 
we shall find that the entire non-trivial dynamical information in the four-loop four-point 
amplitude is contained in a single nonplanar graph; all other graphs are related to this one 
by the duality. 

While a general proof of the duality conjecture is yet to be given beyond tree level, 
the four-loop construction we offer in this paper provides further evidence in favor of it, 
in the form of a highly nontrivial example. In this work, we have confirmed the duality- 
based construction by verifying that the integrand matches a complete (spanning) set of 
generalized unitarity cuts. 

Based on the double-copy structure of supergravity amplitudes, we will give a new rep- 
resentation of the four-loop four-point M = 8 supergravity amplitude. This construction 
provides a direct multiloop confirmation of the double-copy property, because we verify the 
generalized unitarity cuts for the new form of the supergravity amplitude, against the cuts of 
the known expression [10], originally constructed using the KLT relations. We also explore 
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the ultraviolet properties of the amplitude in D > 4 dimensions. An essential feature of the 
new representation is that the UV behavior is manifest: Individual integrals diverge loga- 
rithmically in precisely the same critical dimension D c as their sum. This property did not 
hold for the previous form of the amplitude [10]. The critical dimension is also the same as 
that for the planar and single-trace sectors of M = 4 sYM theory. In a previous paper [10], 
we showed that the supergravity amplitude is finite for D < 11/2, which is also the bound 
obeyed by A* = 4 sYM theory. However, the previous form of the amplitude did not display 
this bound manifestly. To see the cancellation of potential UV divergences, the integrals 
had to be expanded a few orders in powers of the external momenta. The lack of manifest 
UV behavior in that representation made it difficult to carry out the required integration in 
D = 11/2 and to determine whether the amplitude actually does diverge in this dimension. 
With the new form, this task is greatly simplified, allowing us to carry it out here. 

Due to the double-copy construction, the numerators of the integrands for the M = 8 su- 
pergravity amplitude are perfect squares. However, they multiply propagator denominators 
that do not have definite signs. Therefore, individual integrals contributing to the ampli- 
tude can have different signs. To probe whether or not the four-loop amplitude diverges in 
D = 11/2, it is necessary to actually evaluate the UV divergences in the loop integrals in 
this dimension. Using the double-copy form of the four-loop four-point amplitude we do so, 
finding that the M = 8 finiteness bound is in fact saturated at D = 11/2 at four loops, which 
matches the behavior of M = 4 sYM theory. Moreover, we calculate the precise coefficient 
of the A" = 8 supergravity divergence. We find that it exactly matches the coefficient of the 
divergence of the 1 /^-suppressed single-trace term in the four-loop four-point amplitude of 
M = 4 sYM theory, up to an overall rational factor. Although this property is most striking 
at four loops, only emerging after a number of simplifications, it is consistent with lower- 
loop behavior. Presumably this consistent connection is a clue for unraveling the general 
UV properties of N = 8 supergravity. 

Regularization is a crucial point in the construction of loop-level amplitudes in massless 
theories, because such amplitudes are usually either infrared or UV divergent. The issue of 
regularization has been studied in some detail in the context of unitarity cuts in ref. [34], 
where the six-dimensional helicity formalism [35] was suggested as a general means for im- 
plementing either dimensional regularization [36] or a massive infrared regulator equivalent 
to the one in ref. [37]. In the present paper we take advantage of an earlier construction of 
the four-loop four-point amplitude of M = 4 sYM theory, which provides expressions with 
demonstrated validity for D < 6 [11, 34]. In this paper, we compare the D- dimensional uni- 
tarity cuts of the new results with the cuts of the earlier results. We find exact agreement, 
confirming the new representations. 

This paper is organized as follows. In section II we explain our strategy for constructing 
multiloop integrands, illustrating it with the three-loop four-point M = 4 sYM amplitude. In 
section III (and appendix B) we present the new forms of the four-loop integrands of A* = 4 
sYM theory and M = 8 supergravity. We also outline their construction. In section IV, we 
obtain the explicit value of the UV divergence of the N = 8 supergravity amplitude in D = 
11/2 and discuss its properties. We also determine the UV divergence of the color double- 
trace terms in the four-loop A/" = 4 sYM amplitude in D = 6. (We had found earlier [11] 
that the double-trace divergence canceled in the next possible lower dimension, D = 11/2.) 
We give our conclusions and outlook in section V. Several appendices are included. The 
first one gives functional defining relations between the numerators in the four-loop four- 
point amplitude, which are derived from the Jacobi relations after imposing some auxiliary 
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conditions valid for M = 4 sYM theory. Appendix B presents the analytic expressions for 
the numerators. Appendix C gives the values of the vacuum integrals entering the UV 
divergence for the supergravity amplitude in the critical dimension, as well as expressions 
for the vacuum integrals' numerators. Explicit expressions for the color factors for each 
contribution to the full four-loop amplitude may be found online [38], where we also provide 
plain-text, computer- readable versions of the numerator factors and the kinematic dual 
Jacobi relations that they obey. 

II. CONSTRUCTING MULTILOOP INTEGRANDS 

The unitarity method [11, 39-41] has become a general-purpose tool for constructing 
multiloop amplitudes in gauge and gravity theories. In this section we demonstrate how one 
can dramatically reduce the complexity of unitarity-based calculations for gauge theories 
by assuming the conjectured duality between color and kinematics [29, 30]. This duality 
reduces the construction of an amplitude at the integrand level to the determination of 
the numerator factors for a small set of graphs, which we call "master graphs". For the 
four-loop four-point M = 4 sYM amplitude, it suffices to use just two planar master graphs. 
Alternatively, a single nonplanar master graph is sufficient. 

Given the duality-satisfying form of the A/" = 4 sYM amplitude, the M = 8 supergravity 
amplitude can be written down immediately by squaring the A/" = 4 sYM numerator factors. 
We confirm the correctness of the derived gauge and gravity amplitudes using a spanning 
set of generalized unitarity cuts, showing that they agree with our previous forms [10, 11] 
on these cuts. 

A. Duality between color and kinematics 

In general, a massless m-point L-loop gauge-theory amplitude A^m in D space-time 
dimensions, with all particles in the adjoint representation, may be written as 



where g is the gauge coupling constant. The sum runs over the complete set T of m-point L- 
loop graphs with only cubic (trivalent) vertices, including all permutations of external legs. 
In each term, the product in the denominator runs over all propagators of the corresponding 
cubic graph. The integrations are over the independent loop momenta p\. The coefficients 
Cj are the color factors obtained from the gauge-group structure constants by dressing every 
three-vertex in the graph with a factor 



where the hermitian generators of the gauge group are normalized via Tr(T a T fe ) = 5 ab . The 
coefficients rij are kinematic numerator factors depending on momenta, polarization vectors 
and spinors. For supersymmetric amplitudes in an on-shell superspace, the numerators will 
also contain Grassmann parameters. The symmetry factors Si of each graph remove any 
overcount introduced by summing over all permutations of external legs, as well as any 
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FIG. 1: Pictorial Jacobi relation for a group of three graphs. The graphs can represent color 
factors or numerator factors. Except for the connections to the central (pink) lines, the graphs are 
identical in the three indicated by the common (momentum or color) labels a, b, c, d. The 

gray area represents some unspecified subgraph which is identical in all three graphs. 

internal automorphisms — symmetries of a graph with external legs fixed. This form of the 
amplitude may be obtained from representations involving higher-point contact interactions 
(as long as they are built out of / afec, s). One reexpresses all contact terms as the product 
between a propagator and its inverse; i.e. one inserts 1 = p 2 a Jp 2 ail and assigns the inverse 
propagator p^. to be part of a numerator factor n t . The duality conjecture [29, 30] states 
that there should exist a representation of the amplitude where the numerator factors rij 
satisfy equations in one-to-one correspondence with the Jacobi identities of the color factors. 
Explicitly, it requires that 

Ci = Cj + C k => rn = rij + n k , (2.3) 

where the first relation holds, thanks to the usual color Jacobi identity, for any triplet k) 
of graphs which are identical within the gray region in fig. 1. Moreover, the numerator factors 
carry the same antisymmetry properties as color factors, i.e. if a color factor changes sign 
under the interchange of two legs, then so does the corresponding kinematic numerator 
factor: 

d -> -d => m-> -m . (2.4) 

These relations are conjectured to hold to all multiplicities, to all loop orders in a weak- 
coupling expansion, and in any dimension in both supersymmetric and non-supersymmetric 
Yang-Mills theory. Such relations were noticed long ago for four-point tree amplitudes [42]. 
Beyond the four-point level, the relations are rather nontrivial and only work after appro- 
priate rearrangements of the amplitudes. 

While the sign of the complete numerator factor Cjii of a graph is unambiguous, the 
sign of each factor, as well as the signs in eq. (2.3), may be changed as a consequence of the 
relations (2.4) by simply interchanging two lines at any vertex. This ambiguity reflects the 
different possible sign conventions in Jacobi identities. In this paper, the sign of each color 
factor Ci (and implicitly of n*) is fixed by the corresponding graph figure: for each integral, 
Ci is built out of the structure constants f abc corresponding to each trivalent vertex, with 
legs a, b, c ordered clockwise in the plane of the figure. 

The kinematic version of the Jacobi identity (2.3) is the key equation for the duality. At 
loop level, this equation relates graph numerators at the integrand level, as illustrated in 
fig. 1. Therefore it is important to properly line up both external and internal momenta. 



6 



There is one such equation for every propagator of every graph. Of course, many of the 
equations are simply related to each other by automorphic symmetries of graphs, and the fact 
that any given equation can be obtained starting from each of the three contributing graphs. 
Simultaneous consideration of all relations gives a system of linear functional equations that 
the amplitude's numerators should obey. As we will see, only a tiny subset of the possible 
equations needs to be used when solving the system. Once a tentative solution is found, one 
must verify that the full set of equations is satisfied, in order to have a duality-satisfying 
representation of the amplitude. The existence of at least one solution consistent with the 
unitarity cuts is the critical assumption of the conjecture. Indeed, as we will see, the system 
of equations at four loops is quite nontrivial and the emergence of a solution is striking. 

There is by now substantial evidence in favor of the duality, especially at tree level (L = 
0) [43-46] , where explicit representations of the numerators in terms of partial amplitudes 
are known for any number of external legs [47]. A consequence of this duality is the existence 
of nontrivial relations between the color-ordered partial tree amplitudes of gauge theory [29] , 
which have been proven both from field theory [48] and string theory [49] perspectives. These 
relations were important in the recent construction of all open-string tree amplitudes [50]. 
A partial Lagrangian understanding of the duality has also been given [45]. An alternative 
trace-based presentation of the duality relation (2.3), which emphasizes its group-theoretic 
structure, was described recently [51]. 

While less is known at loop level, several nontrivial tests have been carried out. In par- 
ticular, it has been confirmed that the duality holds for the three-loop four-point amplitude 
of Af = 4 sYM theory [30]. (The one- and two-loop four-point amplitudes [8, 52, 53] in this 
theory also manifestly satisfy the duality). Similarly, the duality-satisfying five-point one-, 
two- and three- loop amplitudes of M = 4 sYM theory have recently been constructed [33] . 
The color-kinematic duality is also known to hold [30] for the two-loop four-point identical- 
helicity amplitude of pure Yang-Mills theory [54]. At present there is no proof that the 
system of equations generated by the duality (2.3) always has a solution consistent with the 
unitarity cuts of a given theory, so it needs to be checked case by case. In section III we will 
find a solution for the four-loop four-point amplitude of M = 4 sYM theory. 

Perhaps more surprising than the duality itself is a consequent relation between gauge and 
gravity amplitudes. Once the gauge-theory amplitudes are arranged into a form satisfying 
eq. (2.3), the numerator factors of the corresponding L-loop gravity amplitudes, Aim\ can 
be obtained simply by multiplying together two copies of gauge-theory numerator factors [29, 



where k is the gravitational coupling. The hi represent numerator factors of a second gauge- 
theory amplitude and the sum runs over the same set of graphs as in eq. (2.1). At least 
one family of numerators (rij or hi) must satisfy the duality (2.3). The construction (2.5) 
is expected to hold in a large class of gravity theories, including all theories that are the 
low-energy limits of string theories. At tree level, this double-copy property encodes the 
KLT relations between gravity and gauge-theory amplitudes [31]. For M = 8 supergravity 
both rii and hi are numerators of M = 4 sYM theory. 

The double-copy formula (2.5) has been proven [45] for both pure gravity and for Af = 8 
supergravity tree amplitudes, under the assumption that the duality (2.3) holds in the cor- 
responding gauge theories, pure Yang-Mills and M = 4 sYM theory, respectively. The non- 
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trivial part of the loop-level conjecture is the existence of a representation of gauge-theory 
amplitudes that satisfies the duality constraints. The double-copy property was explicitly 
checked for the three-loop four-point amplitude of M = 8 supergravity in ref. [30] by com- 
paring eq. (2.5) against a spanning set of unitarity cuts of the previously-calculated ampli- 
tude [9, 55]. Here we perform a similar check for the four-loop four-point amplitude of M = 8 
supergravity, using the maximal cut method [41, 55]. The duality and double-copy property 
have also been confirmed in one- and two-loop five-point amplitudes in M = 8 supergrav- 
ity [33]. For less-than- maximal supergravities, the double-copy property has been checked 
explicitly for the one- loop four- and five-graviton amplitudes of N = 4, 5, 6 supergravity [56] 
by showing it matches known results [57]. In the two-loop four-graviton amplitudes in these 
theories, it has been verified to be consistent with the known infrared divergences and other 
properties [58]. The double-copy property also leads to some interesting relations between 
certain M > 4 supergravity and subleading-color sYM amplitudes [59] . 

B. Calculational setup 

We now demonstrate how the conjectured duality between color and kinematics stream- 
lines the construction of integrands of multiloop gauge-theory amplitudes. We first give an 
overview of the procedure and illustrate it with the three-loop four-point amplitude, before 
turning to the four-loop case in the following section. 

To start the construction we enumerate the graphs with only cubic vertices that can 
appear in a particular amplitude. Although this step can be carried out in many different 
ways, we describe one that conveniently also generates the needed duality relations: We 
assume we have a given set of known cubic graphs {e.g. at four loops we can start with the 
planar cubic graphs ones given in ref. [5]). Any missing graphs can then be generated by 
applying the Jacobi relations (2.3) to the set of known graphs. New graphs generated in 
this way are then added to the list of known ones. This process continues recursively, until 
no further graphs or relations are found. At the end of the process all cubic graphs related 
via the duality are known, and all duality relations (2.3) have been written down. 

The next step is to solve the relations thus generated. This is the most complicated part 
of the construction. We can, however, simplify the step by dividing it up into two separate 
parts, the first of which is straightforward. First use a subset of the duality relations to 
express all numerator factors in terms of the numerators of a judiciously chosen small set 
of graphs, which we call "master graphs". We identify master graphs by systematically 
eliminating numerator factors from the duality relations, via a functional analog of the 
standard row reduction of systems of linear equations. This problem is analogous to the 
reduction of loop integrals to a set of master integrals using the Laporta algorithm [60]. In 
both cases, there is freedom to change the order in which the linear equations are solved. 
Here, there is a freedom in the choice of master graphs, which is equivalent to a choice of path 
in solving the system of duality relations. In all cases we have examined, it is convenient to 
choose the master graphs to be planar (although such a restriction does not necessarily yield 
the smallest set). This choice has the advantage that the planar contributions are relatively 
simple and well studied in the literature. In particular, the planar contributions to the 
four- loop four-point amplitude have a fairly simple form [5]. For the three- loop four-point 
M = 4 sYM amplitude we only need a single master graph [30, 32]. In section III, we will 
find that at four loops we can express all numerators in terms of the numerators of only two 
planar master graphs (or a single nonplanar master graph). 
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After the reduction of the system of duality constraints, our task is to find explicit 
expressions for the master numerators. As with any functional equations, a good strategy is 
to write down Ansatze for the master numerators. The Ansatze are then constrained using 
input from unitarity cuts, as well as symmetry requirements on both the master numerators 
and on the numerators derived from them through the duality relations. 

In addition to the duality relations (2.3) and unitarity cuts, we may add extra constraints 
on numerator factors, motivated by our prejudices about the structure and properties of the 
amplitude. Although not necessary, such constraints, when well chosen, can greatly facilitate 
the construction. To find the four-loop four-point M = 4 sYM amplitude we use the following 
auxiliary constraints, which are known to be valid for the duality-satisfying numerators at 
three loops [30]: 

1. One-loop tadpole, bubble and triangle subgraphs do not appear in any graph. 

2. A one-loop n-gon subgraph carries no more than n — 4 powers of loop momentum for 
that loop. 

3. After extracting an overall factor of stA^ ee , the numerators are polynomials in D- 
dimensional Lorentz products of the independent loop and external momenta. 

4. Numerators carry the same relabeling symmetries as the graphs. 

We will assume that these observations carry over to the four-loop four-point amplitude. 
If one of these auxiliary conditions had turned out to be too restrictive, it would have led 
to an inconsistency with either the unitarity cuts or the duality relations. We would then 
have removed conditions one by one until a consistent solution were found. As we shall 
see in section III, these auxiliary constraints are quite helpful for quickly finding a duality- 
satisfying representation for the four-loop four-point amplitude. A surprisingly small subset 
of generalized unitarity cuts is then sufficient to completely determine this amplitude. 

The specific auxiliary constraints that should be imposed depend on the problem at hand. 
The third constraint is clearly specific to the four-point amplitude, and should be modified 
for higher-point amplitudes because they have a more complicated structure. For the five- 
point case, however, a simple generalization has been found [33], involving pre-factors that 
are proportional [44] to linear combinations of five-point tree-amplitudes. For amplitudes in 
less supersymmetric theories, the first and second conditions should be relaxed (in addition 
to the third one), because one- loop triangle and bubble subgraphs are known to appear in 
such theories. 

C. Three-loop warmup 

To illustrate the above procedure in some detail, we reconstruct the well-studied three- 
loop four-point amplitude of M = 4 sYM theory. This amplitude was originally constructed 
in refs. [9, 55]. A form compatible with the duality (2.3) was then found [30]. Here we 
describe how to streamline the construction of the latter form, before following a similar 
procedure at four loops. 

A straightforward enumeration shows that there are 17 distinct cubic graphs with three 
loops and four external legs, which do not have one-loop triangle or bubble (or tadpole) 
subgraphs. Only 12 contribute to the amplitude, as shown in ref. [30]. These 12 nonvanishing 
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FIG. 2: The 12 nonvanishing graphs used in the construction of the J\f = 4 sYM and N = 8 
supergravity three-loop four-point amplitude. The shaded (pink) lines mark the application of the 
duality relation used to determine the numerator of the graph. The external momenta are outgoing 
and the arrows mark the directions of the labeled loop momenta. 

graphs, shown in fig. 2, are sufficient for explaining the construction. (Had we kept all 17 
graphs, the construction would be only slightly more involved, with the result that the 
numerators of the additional five graphs vanish.) 

Each numerator depends on three independent external momenta, labeled by ki,k 2 ,k 3 , 
and on (at most) three independent loop momenta, labeled by / 5 ,/ 6 ,/ 7 , as well as on the 
external states. The Mandelstam variables are s = (&i+&2) 2 ; t = (ki+k^) 2 and u = (ki+k 3 ) 2 . 
We denote the color-ordered tree-level amplitude by A^ ee = A^ ee (l, 2, 3, 4). The four- 
point amplitudes of M = 4 sYM theory are special. Supersymmetry Ward identities fix the 
external state dependence, and imply that an overall prefactor of A^ ee , or equivalently the 
crossing-symmetric prefactor stA^ ec [9, 55], can be extracted from every numerator factor 
rS x \ leaving behind new numerator factors that depend only on the momenta, 

n {x) = sL4* ree (l,2,3,4)iV^, 

N ( X ) = N {x) (k 1 ,k 2 ,k 3 ,l 5 ,l 6 ,l 7 ) . (2.6) 

Here (x) refers to the label for each graph in fig. 2. This result has been argued to be valid 
in any dimension D < 10 [11], justifying the third assumption above for these amplitudes. 
The homogeneity of the Jacobi relations implies that they hold for N^ x ' just as for rS x \ The 
crossing symmetry of stA^ ee implies that the symmetry properties of are the same as 
those of rS x \ 
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Next, we write down a subset of the duality relations that allows us to identify the 
master graphs [32]. For the three- loop four-point amplitude, one such restricted set of 
duality relations is: 

jy(a) = N ^( kl ,k 2 ,h,k,l 6 ,l 7 ), 

N (h) = N w(k u k 2 ,k 3 ,hJ 6 J 7 ), 

N {c) = AT (a) (fci,A; 2 ,fc 3 ,M 6 ,Z 7 ), 

N (d) = N W( k3) kl} fc 2j l 7) k) fcli3 + l 7 ) + N^(k 3 , k 2 , h, l 7 , l 6 , fc 2j3 + h- h) , 

iv(f) = N^ik^hMMM , 

Arte) = N^ik^kMM, 

AT(h) = -N( s \ki, k 2 , k 3 , h, l 6 , fc 1)2 -k- l 7 ) - N^(h, k 3 , k 2 , k - k, h-k + h- h, 2 , k) , 

N {i) = N (°)( kl , k 2 , k 3 , h, l 7 , l 6 ) - N®(k 3 , k 2 , h, -h -k- k, -h - h, h) , 

AtO) = iV^(fci, k 2 , k 3 , / 5 , / 6 , l 7 ) - A r(e) (/c 2 , fci, A; 3 , / 5 , / 6 , / 7 ) , 

atW = N^^krMMMMM - NV\k 2 MMMM,h) , 

N (D = Nte(h,k 2 ,k 3 ,k,lM - N^(k 2 ,h,k 3 ,l 5 ,l 6 ,l 7 ) , (2.7) 

where fcjj = fcj + fcj. For convenience we have suppressed the canonical arguments 
(ki, k 2 , k 3 , 15, hi h) of the numerators on the left-hand side of the equations (2.7), as we 
will often do in the remainder of the paper. Each relation specifying an is generated by 
considering the dual Jacobi relations focusing around the lightly colored (pink) line labeled 
J x in fig. 2. In general, the duality condition relates triplets of numerators; sometimes, 
however, one or two of the numerators vanish because the associated graph has a one-loop 
triangle subgraph forbidden by our auxiliary constraints. Specifically, for five of the above 
equations, the duality sets pairs of numerators equal; this occurs because the third term in 
the triplet of numerators of eq. (2.3) vanishes due to the presence of a triangle subgraph. 
The above system can be used to express any numerator factor in terms of combinations of 
the numerator N^ c \ with various different arguments. Thus, graph (e) can be taken as the 
sole master graph. This is a convenient choice, but not the only possible one; for example, 
either graph (f) or (g) can also be used as a single master graph. None of the remaining 
nine graphs, however, can act alone as a master graph. 

One valid numerator factor (consistent with unitarity cuts) for graph (e) is the "rung-rule" 
numerator [53], 

= s {h + k 4 f. (2.8) 

With this numerator, the graph possesses dual conformal symmetry. However, it turns out 
that this numerator is incompatible with the duality between color and kinematics (2.3). 

We are therefore looking for a modification of consistent with both the maximal cut 
of the graph and with the duality constraints (2.7). We start by requiring that the maximal 
cut of graph (e) is correct, and that the auxiliary constraints in section II B are satisfied. 
That is, the numerator iV( e ) has mass dimension four and possesses the symmetry of the 
graph; no loop momentum for any box subgraph in (e) appears in it (ruling out 1$ and l 7 ); 
and is at most quadratic in the pentagon loop momenta Z5. (This last condition is 
looser than the second auxiliary constraint in section II B; we will tighten it shortly.) The 
symmetry condition implies that iV( e ) should be invariant under 

{ki k 2 ,k 3 h,l 5 fci + k 2 - l 5 } . (2.9) 



11 



The most general polynomial consistent with these constraints is of the form, 

A (e) = s(k + h) 2 + (as + pt)l\ + (js + St)(k - h) 2 + (as + (3t)(k - h - k 2 ) 2 , (2.10) 

where the four parameters a, /3, 7, 5 are to be determined by further constraints. All added 
terms are proportional to inverse propagators and therefore vanish on the maximal cut. 
Thus, since eq. (2.8) is consistent with the maximal cuts so is eq. (2.10). 

According to the second auxiliary constraint in section II B, the numerator of a pentagon 
subgraph should be at most linear in the corresponding loop momentum, not quadratic 
as assumed above. Therefore the coefficient of I 2 in eq. (2.10) should vanish, yielding the 
relations 7 = —1 — 2a and § = —2/3. Consequently, the Ansatz for is reduced to 

iV (e) = s(t 45 + r 15 ) + (as + (3t)(s + r 15 - r 25 ) , (2.11) 

where we use the notation, 

Tij = 2k t ■ lj (i < 4, j > 5) , Tij = 2h ■ I, (i, j > 5) . (2.12) 

Now there are just two undetermined parameters, a and /3. 

To determine one of the remaining parameters we use the properties of graph (j), and 
the expression for its numerator in terms of the numerator of graph (e), which is given by 
the 9th duality constraint in eq. (2.7). Inserting eq. (2.11) into this relation leads to 

A^ j ) = s(l + 2a- 0)(t 15 - r 25 ) + f3s(t - u) . (2.13) 

Because the smallest loop in graph (j) carrying l 5 is a box subgraph, our auxiliary constraints 
require that this momentum be absent from N^ e \ Setting the first term in eq. (2.13) to zero 
implies that /3 = 1 + 2a, which in turn leads to 

= s(t A5 + 7i 5 ) + (a(t -u)+ t)(s + r 15 - r 25 ) , (2.14) 
A (j) = (l + 2a)(t-u)s, (2.15) 

leaving undetermined a single parameter a. 

There are a variety of ways to determine the final parameter. For example, one can use 
planar cuts to enforce that the planar part of the amplitude is correctly reproduced. A 
particularly instructive method is to use the duality relations to obtain the numerator for 
planar graph (a) in terms of master numerator N^ c \ The numerator is quite simple 
once we impose the condition that a one-loop box subgraph cannot carry loop momentum. 
Since three independent one-loop subgraphs of graph (a) are boxes, the numerator 
cannot depend on any loop momenta. Indeed, the iterated two-particle cuts, or equivalently 
the rung insertion rule [53], immediately fix this contribution to be 

A (a) = s 2 . (2.16) 

Solving the duality relations (2.7) to express in terms of we find, 

jV(a) = N (-)( ku k 2 ,k 4 ,-k 3 + k-l 6 + l 7 ,l 5 -l 6 ,-l 5 ) 
+ N {c \k 2 , k ly fc 4 , -h -k + h, ~h, h - k) 

- A (c) (fc 4j k u k 2 , l 6 - l 7 , l 6 , l 5 - l 6 ) - A (c) (fc 4 , k 2 , h,l 6 - l 7 , / 6 , -l s ) 

- (k 3 , h, k 2 , h, k, h - k) - A (e) (k 3 , k 2 , h, l 7 , k, ~k) ■ (2.17) 
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Plugging in the value of the numerator factor in eq. (2.14), we obtain 

= s 2 + (l + 3a) ((Tie - r m )s - 2(r 17 + r 37 )s + (r 16 - 2r 17 - r 26 + 2r 27 )t + 4ut) . (2.18) 

Demanding that this expression matches the numerator factor given in eq. (2.16), 
or alternatively that it is independent of loop momenta, fixes the final parameter to be 
a = —1/3. This constraint completely determines the numerator of graph (e) to be 

= S (r 45 + r 15 ) + -{t - s){s + r 15 - r 25 ) , (2.19) 

matching the result of ref. [30]. 

Remarkably, numerator (e) in eq. (2.19) generates all other numerators N^ x \ via eq. (2.7), 
giving us the entire integrand at three loops. For all graphs, the resulting numerators re- 
produce the expressions quoted in ref. [30], and the resulting amplitude matches previous 
expressions [9, 55] on all D-dimensional unitarity cuts. As already noted, it is highly non- 
trivial to have a consistent solution where all duality relations hold, all numerators have the 
graph symmetries and all unitarity cuts are correct. Squaring these numerators N^ x \ using 
eq. (2.5), immediately yields the numerators for the three-loop four-point M = 8 supergrav- 
ity amplitude. This form has also been confirmed against previous expressions [9, 55] on a 
spanning set of D- dimensional unitarity cuts [30]. 

We shall use the same streamlined strategy to construct the four-loop four-point ampli- 
tude in section III. Before carrying out this construction, however, we need to address an 
important subtlety that appears in the construction of the three-loop amplitude and affects 
the four-loop construction as well. 

D. Comment on one-particle-reducible graphs and snails 

Beyond tree-level, the on-shell three-point amplitudes of M = 4 sYM theory vanish. The 
appearance of one-particle reducible (1PR) graphs in the three- loop four-point amplitude 
may therefore seem surprising. Indeed, graphs (i), (j) and (k) of fig. 2 do not appear in the 
original representations of the same amplitude [9, 55] . The existence of 1PR graphs may seem 
to imply that the three-point amplitude is non-vanishing. However, these graphs' numerators 
are proportional to the Mandelstam invariant s, which is also the inverse propagator for the 
sole line on which the graph is 1PR. Thus, the superficially 1PR graphs are in fact just one- 
particle-irreducible (1PI) contact graphs. Even though they are kinematically equivalent to 
1PI graphs, the non-contact form of graphs (i), (j) and (k) in fig. 2 is needed to describe easily 
their color structure, and to allow the amplitude to obey the duality (2.3) between color 
and kinematics. As we shall see, this feature continues at four loops, where we encounter, 
not only graphs with three-point subgraphs, but also nontrivial two-point subgraphs. Some 
of these graphs contain four-loop two-point bubble subgraphs on external legs, and must be 
treated with particular care. 

At first sight, it may appear surprising that two- and three-point subgraphs show up; 
indeed in M = 4 sYM theory we expect the vanishing of on-shell two- and three-point loop 
amplitudes. This property has been known in string theory for some time [61]. By taking 
the low-energy limit, it should hold in field theory as well. A direct field theory argument 
for the vanishing of the on-shell two-point function can be made as follows: Quite generally, 
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a (diagonal) on-shell two-point loop contribution represents a correction to the mass of the 
corresponding field. Gauge invariance forbids such a term from being generated in the gluon 
two-point function. The chirality of M = 4 sYM gluinos forbids such mass terms from being 
generated by perturbative quantum effects for fermions as well. Thus, gluon and gluino 
two-point functions vanish on shell. Manifest off-shell M = 1 supersymmetry, which can be 
maintained, then implies that the scalar field two-point function also vanishes on shell. 

We can also argue that three-point amplitudes vanish on shell. Because M = 4 super- 
symmetry relates all such amplitudes to each other, it suffices to focus on the scattering 
amplitude of two fermions and one scalar. Up to an SU (4) i?-symmetry transformation, we 
may further assume that neither of the fermions is the M = 1 superpartner of the gluon. 
Thus we consider only the interaction between M = 1 matter multiplets. Conservation 
of the matter i?-symmetry subgroup 577(3) requires that the three-field interaction is con- 
trolled by S77(3) invariance, and thus is either holomorphic or antiholomorphic. Now, in 
the effective action language, the three-point amplitude originates either from terms in the 
superpotential or the Kahler potential. Due to the perturbative nonrenormalization of the 
superpotential [62], only the tree amplitude comes from the former. A nonvanishing loop 
amplitude can only originate from a correction to the Kahler potential. For this case, a 
nonvanishing full superspace integral and Lorentz invariance require that the product of 
three chiral superfields containing the relevant wave functions must be accompanied by at 
least two additional superderivatives. In turn, this implies that the product of one scalar 
and two fermion wave functions is always accompanied by an external momentum invariant, 
originating from the superspace integration measure. For massless fields, any such product 
vanishes on shell. Thus, all quantum corrections to three-point amplitudes in M = 4 sYM 
theory vanish on shell, completing the argument. 

While these arguments confirm the vanishing of two- and three-point amplitudes at one- 
loop and beyond, we emphasize that this does not mean that we cannot have graphs with 
two- and three-point subgraphs. However, when such graphs appear they should always 
carry factors that make their contributions vanish whenever legs are cut (placed on shell) 
to isolate two- and three- point subamplitudes. Indeed, we shall find that at four points, 
through four loops, all such graphs with three- or four-point subgraphs can be absorbed 
as contact terms in other graphs. This property is consistent with the fact that previous 
representations of the three- and four- loop amplitudes [9, 11, 55] do not use any 1PR graphs 
with two- or three-point loop subgraphs. 

At L loops, the four-point amplitude in M = 4 sYM theory is expected to have a repre- 
sentation with at most 2(L — 2) powers of the loop momentum in the numerator of each 1PI 
cubic graph [8, 63]. At three loops, this power counting allows for cancellation of one inter- 
nal propagator, as in graphs (i), (j) and (k) of fig. 2. However, it precludes the existence 
of two-point graphs or propagator corrections (and tadpole graphs), which would require 
two inverse propagators or four powers of the loop momentum in the numerator. On the 
other hand, at four loops and beyond, such graphs can and indeed do appear. Propagator 
corrections can be of two types: 

1. on internal legs, as shown in fig. 3(a), or 

2. on external legs, as shown in fig. 3(c). 

In both cases the graph's numerators must contain momentum invariants that cancel out 
the unwanted poles, so that they are kinematically equivalent to the 1PI graphs shown in 
fig. 3(b) and fig. 3(d), respectively. 
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FIG. 3: Graphs (a) and (c) are generic propagator correction graphs that can appear at four 
loops and beyond if we have a cubic organization of graphs. Graphs (b) and (d) are rewritings of 
these graphs, which make explicit that in M = 4 sYM theory numerator factors always cancel the 
propagators that are external to the loops in the four-point amplitude. 

For case 1, this cancellation is straightforward because the momentum invariant is non- 
vanishing for generic on-shell kinematics. For case 2, the external leg corrections, the mech- 
anism is more subtle. On the one hand, because the amplitudes have on-shell external legs, 
a propagator in fig. 3(c) diverges: l/(ki + k% + k^) 2 = 1/0. On the other hand, from the 
vanishing of the on-shell two-point function we expect that the numerator of fig. 3(c) is 
proportional to k\ = 0. In order to resolve this 0/0 ambiguity, we need to regulate the 
external leg by taking k\ ^ 0, and cancel factors of k\ between numerator and denominator. 
This procedure yields the "snail graph" 1 in fig. 3(d), which is perfectly well behaved at the 
level of the integrand, even with all external momenta on-shell. 

It is important to note that the snail graph in fig. 3(d) contains a scale- free integral, 
which vanishes by the usual rules of dimensional regularization. We cannot, however, simply 
ignore these contributions. In dimensional regularization, scale-free integrals evaluate to zero 
because of cancellations between infrared and UV singularities. Ignoring the snail graph 
contributions in M = 4 sYM theory would lead to incorrect values for the UV divergences. 2 
Since we are interested in this paper in the coefficient of the UV divergences, these snail 
graphs must be included. 

While enforcing the duality constraints (2.3) brings the phenomenon of snail graphs to the 
forefront, we emphasize that the potential appearance of such contributions to amplitudes 
is independent of the color-kinematic duality. Snail contributions can in principle occur 
within any representation; it is therefore important to always check the unitarity cuts for 
such contributions. Because snails are associated with external leg contributions, ordinary 
unitarity cuts fail to detect them, and generalized cuts are required. The momenta of the 
states crossing the cut must either be complex, or else have an indefinite sign of their energy. 



1 With suitable imagination, the graph resembles a snail (as much as a penguin diagram resembles a 
penguin). 

2 In QCD, propagator corrections on external legs can be ignored because the UV divergences are known 
a priori. It is therefore quite simple to restore the missing terms. In J\[ = 4 sYM theory, UV divergences 
in D > 4 are unknown a priori. 
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(We have examined such cuts, and have confirmed thereby that no snail contributions are 
present in the representation found in ref. [11].) 

Although the snail contributions are important for M = 4 sYM amplitudes, they will not 
infect the corresponding Af = 8 supergravity amplitudes. This may be understood heuristi- 
cally as a consequence of the double-copy formula (2.5). In Af = 4 sYM theory, graphs of 
the form in fig. 3(c) carry a factor of in their numerator to cancel the 1/0 from the propa- 
gator. In Af = 8 supergravity we get a second factor of from the second copy, making the 
numerator vanish faster than the denominator, and giving a vanishing snail contribution. 
Below we confirm this heuristic argument directly from unitarity cuts. 

Finally, we remark that very similar considerations appear in the analysis of inverse 
derivative factors arising from the collision of vertex operators in the discussion of nonrenor- 
malization conditions for amplitudes in superstring theory — see section 3.2 of ref. [64]. 



III. THE FOUR-LOOP FOUR-POINT INTEGRAND 

We now turn to the construction of the four- loop four-point amplitude and follow the same 
strategy as described in the previous section for the corresponding three-loop amplitude. 



A. Overview of the result 

We will find that, in terms of the 85 distinct graphs of figs. 5-11, the four-loop sYM 
amplitude is given by 

Af = «« e t /(n (gfe) iWf 1 • (3-D 

S 4 i=l J \j'=5 V > 7 * VYa i =\Pcx l 

where l 5 ,l 6 ,l 7 , l s are the four independent loop momenta and ki,k 2 , k 3 are the three indepen- 
dent external momenta. The p ai are the momenta of the internal propagators (corresponding 
to the internal lines of each graph i), and are linear combinations of the independent loop mo- 
menta lj and the external momenta k m . In the case of 1PR graphs, some p ai will depend only 
on the external momenta. As usual, d D lj/(2ir) D is the D-dimensional integration measure 
for the j th loop momentum. The numerator factors N^kj, lj) are polynomial in both internal 
and external momenta, and are given in appendix B. The color factors Cj = Cf ia2a3a " 1 are 
collected online [38], but they can also be read directly off the figures. The full amplitude 
is obtained by summing over the group £4 of 24 permutations of the external leg labels. 
Overcounts are removed by the symmetry factors Si, which include both external symmetry 
factors (related to the overcount from the sum over <S 4 ), as well as any internal symmetry 
factors associated with automorphisms of the graphs holding the external legs fixed. As at 
three loops, we extract the crossing-symmetric, ^-invariant prefactor stA^, which con- 
tains all dependence on the external states. (Notice that we have used a slightly different 
notation for the independent loop momenta lj in eq. (3.1), compared with pi in eq. (2.1).) 

Out of the 85 integrals in eq. (3.1), graphs 50 and 79 are somewhat peculiar: Their 
integrands are nonvanishing, but they integrate to zero. The vanishing of their integrals 
can be seen from symmetry considerations alone. For example, graph 50 has a symmetry 
exchanging legs 1 and 4, and legs 2 and 3, flipping the graph across a vertical midline. It is 
easy to check that the color graph C50 picks up a minus sign under this operation; therefore 
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the kinematic integrand must also be antisymmetric, causing the integral to vanish. In fact, 
as the duality between color and kinematics might suggest, the color factors C50 and C79 
vanish after the internal color sum is carried out (for any gauge group G). However, both 
graphs give nonvanishing contributions to the M = 8 supergravity amplitude. Therefore 
we retain them here. (While the vanishing gauge-theory integrals are odd under the above 
relabeling of the loop momenta, the double-copy property makes the gravity integrals even 
under the same relabeling.) 

As we discussed in section II D, graphs 83-85, (displayed in fig. 11), superficially appear 
as propagator corrections on external legs. These graphs give rise to the snail contributions 
described there, after an external propagator is canceled by a corresponding factor in the 
numerator. 

Using the double-copy relation (2.5), the four-loop four-point M = 8 supergravity am- 
plitude is obtained simply by trading the color factor for hi = sL4 4 ree iVj in eq. (3.1). 
Employing the relation s 2 t 2 A£ ce A£ ee = istu M"£ ee and changing the gauge coupling to the 
gravitational coupling, we have 

-?Hi)'VEi:/(ni)^, <w> 

Sa i=l V i=5 v ' 7 l i.lai=lfai 

where Ni(kj, lj) are the gauge-theory numerator factors given in appendix B. In contrast to 
the sYM amplitudes, potential snail contributions from graphs 83-85 vanish identically, as 
expected from our heuristic argument in section II D, and confirmed by an analysis of the 
unitarity cuts. 

The symmetry factors appearing in eqs. (3.1) and (3.2) are given explicitly as: 

85 . 

TT-fi = jh + jh + + jh + + \h + \h + -^8 + \h + \ho + \hi + \l\2 + \hz 
i=i ^ 

+ |/l4 + 4/15 + Il6 + 2^!7 + ^18 + Il9 + -^20 + hi + I22 + I23 + |-^24 + ^25 + ^26 



+ 2^27 + 4-^28 + ^29 + |-^30 + \hl + ^32 + ^33 + |-^34 + ^35 + ^36 + f-^37 + |-^38 
+ 2^39 + 4-^40 + |-^41 + -^42 + I 43 + |-^44 + 4X45 + §-^46 + |-^47 + |-^48 + |-^49 + f-^50 
+ 2^51 + ^52 + 4-^53 + 2^54 + 2^ 55 + f-^56 + f ^57 + |-^58 + f-^59 + |-^60 + 4-^61 + |-^62 
+ |-^63 + 2^64 + f-^65 + \ -^66 + |-^67 + f-^68 + |-^69 + |-^70 + f-^71 + §-^72 + |-^73 + 4 -^74 
+ 4-^75 + 2^76 + f-^77 + 4" -^78 + g-^79 + 16^80 + f-^81 + 7^82 + 4-^83 + |-^84 + ^85 , 

(3.3) 

where Ii should be interpreted only as placeholders for the graphs, including both the nu- 
merator or color factors, in either theory. 



B. The calculation 

As at three loops, the construction of the amplitude begins by writing down a sufficient 
number of duality constraints so that a set of master numerators can be identified. We 
have constructed a set of duality equations similar to the three-loop ones of eq. (2.7). In 
appendix A we collect a set of simplified equations derived from these duality constraints 
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(18) (28) 

FIG. 4: The planar master graphs, 18 and 28. The numerators and color factors of all other graphs 
are generated from the numerators and color factors of these two graphs through kinematic Jacobi 
relations. 



by imposing the four auxiliary constraints presented in section II B. Because of these addi- 
tional simplifications, the duality equations in the appendix are valid only for M = 4 sYM 
theory. The duality equations allow us to express all non-snail numerators directly as linear 
combinations of the numerators iV 18 and iV 28 . The corresponding graphs are shown in fig. 4; 
we will choose them as the master graphs. 

It is interesting to note that as an alternative we can use a single nonplanar master 
graph that does the same job, such as graph 33 (or equivalently 35 or 36, which have 
identical numerators up to a sign). However, we prefer to use planar graphs as master graphs 
because their numerators have a somewhat simpler structure. If we choose planar graphs 
as master graphs then the minimal number is two. In our treatment the snail contributions 
are only given partially in terms of the master numerators, because the latter are specified 
using on-shell external kinematics, whereas the numerators of the former require an off-shell 
regularization to be nonvanishing. 

Our next task is to determine the master numerators. To this end we begin by construct- 
ing an Ansatz for the numerator factors iVig and N28 that satisfies the auxiliary constraints 
discussed in section II B and a restricted set of duality relations. We then constrain the 
Ansatz by demanding that other duality relations are satisfied, and that the numerator fac- 
tors of the other integrals obey the auxiliary constraints. For both graphs, the numerator 
must be independent of loop momenta I7 and /§ because they are assigned to one-loop box 
subgraphs, whose momenta should not appear in their numerators. Similarly, momenta 1$ 
and l 6 are assigned to one-loop pentagon subgraphs, so N 1S and N 28 should be no more 
than linear in these momenta, according to our auxiliary constraints. Thus, each of the two 
master numerators should be a polynomial built from the monomials, 

M = {s 3 , St 2 , S 2 t, t 3 , Ti^S 2 , Ti^t 2 , T i5 st, T i6 S 2 , T^t 2 , TigSt, TffiTjQS, Ti^Tj 6 t, T§qS 2 , T 56 t 2 , T^st} , 

(3.4) 

where % — 1, 2, 3 labels the three independent external momenta, k±, k 2 , k^. In total this gives 
us a polynomial with 43 terms for each master graph. Labeling the monomials consecutively 
as Mj, and including arbitrary coefficients, we have as our starting Ansatz, 

43 43 

N 18 = J2 a i M i . N 28 = J2 b J M i ■ ( 3 - 5 ) 
3=1 3=1 

The 86 free coefficients Oj and bj are to be determined from various consistency conditions 
obtained from the color-kinematic duality, graph symmetries and unitarity cuts. The number 
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(1) (2) (3) 




FIG. 5: Cubic graphs 1 to 11 that contribute to the four- loop four-point amplitude of N = 4 sYM 
theory and N = 8 supergravity. The labels 1 to 4, indicate the legs carrying external momenta k\ 
to &4. The labels 5 to 8 indicate the propagators carrying the independent loop momenta 1$ to Is- 
The arrows indicate the direction of the momentum. The graphs also specify the color factor of 
the graph, simply by dressing each cubic vertex with an f abc , respecting the clockwise ordering of 
lines at each vertex. 

of free parameters that need to be determined in this construction is remarkably small, 
considering the expected analytic complexity of amplitudes at four loops. 

Using the 86-parameter Ansatz and the solution to the restricted set of duality constraints 
listed in appendix A gives us expressions for the numerator factors of any of the 82 non- 
snail graphs appearing in the amplitude. (The snail graphs will be determined below in 
terms of the non-snail graphs using generalized unitarity cuts.) These expressions do not 
yet satisfy all duality constraints; thus far we have imposed only the relatively few relations 
in appendix A sufficient to determine all numerators in terms of the master numerators, but 
we have not yet accounted for the complete set of duality relations. To further constrain the 
master Ansatz we could require that all other dual Jacobi relations are satisfied; there are on 
the order of 13 x 85 such functional relations (not all independent). An alternate strategy, 
which we follow here, is to first impose the consistency constraints on the numerators of the 
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(12) (13) (14) 




(15) (16) 




FIG. 6: Cubic graphs 12 to 23 that contribute to the four-loop four-point amplitude of N = 4 sYM 
theory and N = 8 supergravity. The labeling is the same as for fig. 5. 

graphs derived from N\% and A^8 through eqs. (A4) and (A5). After obtaining a complete 
solution for all 86 parameters appearing in the ansatz, we then verify that they indeed satisfy 
all remaining duality relations and unitarity cuts. An advantage of this strategy is that it 
allows us to illustrate the remarkably small number of unitarity cuts needed to find the 
complete amplitude, including nonplanar contributions. 

As we shall see, to construct the complete amplitude we need only information about the 
unitarity cuts of the four-loop planar amplitude, obtained previously in refs. [5, 11]. The list 
of constraints needed to fix all 86 parameters in the Ansatz, thus determining the amplitude, 
is remarkably short. It is sufficient to enforce: 

1. the graph automorphism symmetries on numerators N12, and N28', 

2. the maximal cut of graph 12; 

3. the next-to-maximal cut of graph 14, where I5 is the off-shell leg. Graph 68 also 
contributes to this cut. 
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Strikingly, only two rather simple planar cuts are needed to fully determine the amplitude. 
Let us now discuss some details of fixing the parameters. 

We start by analyzing the consequences of the symmetries of the master graph 28: This 
graph is invariant under two independent transformations: 

{h o k 3 , l 5 k 2 - h, h->ki- k, h -H- / 8 } , (3.6) 
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(39) (40) (41) 




(48) 4 1 (49) (50) 



FIG. 8: Cubic graphs 39 to 50 that contribute to the four-loop four-point amplitude of N = 4 sYM 
theory and N = 8 supergravity. The labeling is the same as for fig. 5. 

and 

{k 2 <-> h, h <-> k, h-th- i 7 , h^h- 1 8 } . (3.7) 

Imposing the invariance of numerator N 28 under eqs. (3.6) and (3.7) reduces the number of 
its unknown coefficients from 43 to 14. The other master graph, graph 18, does not have any 
such automorphism relations; we are therefore left to determine a total of 57 parameters. 

We then impose similar symmetry conditions on N\ 2 , which may be written in terms of 
Ni 8 and N 28 as 

N 12 = -N 18 {k 4 , /c3,/ci,/ 6 ,-/ 5 ,-/ 6 ,/ 8 ) + N 18 (k 4 , k 3 , k 2 ,l 6 , k 2 + l 8 ,hJi) 

+ N 18 (k 4 , k 3 , k 2 , k 3 + l 8 , h ~ k, ~h, k) ~ N 28 (k u k 2 , k 3 , l 5 - l 8 , k 3 - l 6 + l 8 , -l 6 , l 8 ) 

+ N 28 (k 2 , k±, k 3 , — Z5, 0, — Iq, l 8 ) — N 28 (k4, k 3 , k%, Iq — l 8 , k 2 — l§ + l 8 , l 8 , l 8 ) 

+ N 28 (ki, k 3 , k 2 , k 3 , ki + l 5 , -k 3 + l 6 , l 8 ) , (3.8) 

by combining the 2nd, 6th, 14th and 21st relations in eq. (A4) in appendix A. Invariance 
under the automorphisms of graph 12 fixes 37 parameters, leaving undetermined 20 param- 
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FIG. 9: Cubic graphs 51 to 65 that contribute to the four-loop four-point amplitude of N = 4 sYM 
theory and M = 8 supergravity. The labeling is the same as for fig. 5. 



eters. Similarly, imposing the graph symmetry condition on N u reduces the total number 
of unknown parameters to 17. (These parameter counts are for the specific set of duality 
relations given in eq. (A4). Using another set of relations would result in somewhat different 
parameter counts; however, the final solution would be the same.) 

We could continue imposing more symmetry constraints on other numerators, but we 
already have a very small set of undetermined parameters. Ultimately, dynamical informa- 
tion provided by unitarity cuts should become necessary. Therefore we will now inspect 
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FIG. 10: Cubic graphs 66 to 79 that contribute to the four-loop four-point amplitude of M = 4 
sYM theory and W = 8 supergravity. The labeling is the same as for fig. 5. 

some cuts. A good starting point is the maximal cut of graph 12. Its explicit value is easily 
obtained using the simple rung-rule numerator of that graph [5, 53], 



N? 2 = s (h + k + h + h) 2 . 



(3-9) 
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FIG. 11: Cubic graphs 80 to 85 that contribute to the four-loop four-point amplitude of M = 4 
sYM theory. Graphs 83, 84, and 85 vanish identically for N = 8 supergravity, but carry an UV 
singularity in the D = 11/2AA = 4 sYM case. The labeling is the same as for fig. 5. 

The rung rule was originally designed to reproduce iterated two-particle cuts. Since maximal 
cuts can be obtained from iterated two-particle cuts by imposing additional cut conditions, 
the rung rule reproduces the maximal cuts as well. Our task is to match N 12 , as obtained 
from the duality relations, and iYfJj in eq. (3.9) on the maximal cut kinematics that uniquely 
single out this graph, i.e. we impose pf = on all 13 propagators of graph 12. Solving these 
conditions, we obtain 



k = , Ti5 = r 25 = r 36 = 

r 38 = r 68 j 728 = r 58 j T 27 : 



Tl6 = -T 26 , Tir = -T 57 , 

-t 37 + r 57 + r 67 , r 18 = — s + r 37 



r 58 — T 67 



(3.10) 

r 78 • 



Thus, on the maximal cut, N" 2 becomes 



A rmax. cut 
iV 12 



s 2 (t - r 26 - r 35 + r 56 ) 



(3-11) 



max. cut 
12 



Requiring that the numerator obtained from the Ansatz (3.5) via eq. (3.8) matches N'- 
on the maximal cut, the number of undetermined parameters is reduced from 17 to 8. 

Finally, all remaining parameters can be determined by requiring the next-to-maximal 
cut of graph 14, where all propagators except for 1/Z§ are placed on shell, to be satisfied. 
Graph 68 also contributes to this cut since it contains the same set of cut propagators. 
Relabeling graph 68 so it matches graph 14, and appropriately weighting the numerators by 
the remaining off-shell propagators, we find under the cut kinematics, 



s (k + h + k 3 y 



(3.12) 



The right-hand side of this equation is the numerator of graph 14, as constructed using 
the rung rule; in the rung-rule representation of the planar four-loop amplitude, graph 68 
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FIG. 12: The snail contributions in a form resolving the 0/0 ambiguity of graphs 83, 84 and 85. 




FIG. 13: A cut that determines the snail contribution in fig. 3(d). In N = 4 sYM theory this cut 
vanishes at the level of the integrand. 

vanishes [5]. This requirement fixes all remaining eight coefficients, giving us a unique 
expression for the master numerators, 

iV 18 = i(6M 2 r 25 +M(2s(5r 25 + 2r 26 )-r 15 (7r 16 + 6t)) 

+ £(ti5T 2 6 - t- 25 (ti 6 + 7r 26 )) + s(4r 15 (t - r 26 ) + 6t 36 (t 35 - r 45 ) 

- n 6 (4t + 5r 25 ) - t 46 (5t 3 5 + r 45 )) + 2s 2 (t + r 26 - r 35 + r 36 + r 56 )) , (3.13) 
^28 = |(s(2ri 5 i + ri 6 (2t - 5r 25 + r 35 ) + 5r 35 (r 26 + r 36 ) + 2t(2r 46 - r 56 ) - 10wt 25 ) 

- 4s 2 r 25 - 6n(r 46 (t - r 25 + r 45 ) + r 25 r 26 ) - t(ri 5 (4r 36 + 5r 46 ) + 5r 25 r 36 )) . 

It is quite striking that the solution we obtain is unique and relatively simple. There are 
presumably other solutions to the duality relations for this amplitude, but finding them 
would require relaxing some of the auxiliary constraints described in section II B. 

This construction completely fixes the values of all numerators from N\ to A*s 2 (subject 
to on-shell external kinematics) using the duality relations in appendix A. 

C. Resolving the snails 

In the previous subsection we obtained the numerators for the master graphs, which 
determine all numerator factors through duality relations. However, it does not resolve 
the 0/0 ambiguity appearing in graphs 83-85 in fig. 11, which needs at least one external 
momentum off shell to be properly defined, if we insist on representing the result in terms of 
graphs with only cubic vertices. As explained in section II D, these contributions are indeed 
finite, as can easily be confirmed explicitly using generalized unitarity cuts, as we do below. 

First we consider the planar snail contributions to the four-loop four-point amplitude. 
To determine them, we evaluate the unitarity cut shown in fig. 13. In the planar limit, the 
integrals contributing to this cut are shown in fig. 14. Because the cut contains an on-shell 
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FIG. 14: The graphs contributing to the cut in fig. 13. The sum of these contributions to the cut 
vanishes in N = 4 sYM theory. 



three-loop three-point subamplitude it must vanish, following the arguments in section II D. 
Therefore we can determine the contribution of the snail graph 83 in terms of that of graphs 
51 and 55. Dressing the numerators of the first three graphs in fig. 14 with the appropriate 
ratios of propagators, we find that 



(L — k] — ko) 2 . 

^51 + — — UV 55 



t 



N, 



55 



fci-H-A;3,Z5— >fci+fc2— 



27 



{h-hfih-h-hfu 



(3.14) 



where we used the relations l\ — > and T45 — > 0, valid on the cut. Here we relabeled the 
momenta of the two contributions of graph 55 (compare to fig. 9) to match the labels of 
graph 51. The two factors on the right-hand side of eq. (3.14) that depend on Z5 are simply 
inverse propagators belonging to graph 51. In order to compare to graph 83 we should 
remove these factors, because it does not have these propagators. Thus, in order for the 
sum of contributions in fig. 14 to vanish, we must have 



-AW 

s 



-( 



fcl-flfc3 



27 



-u 



(3.15) 



where we start with the labeling in fig. 12 and relabel accordingly. (Because the snail nu- 
merators turn out to have no loop-momentum dependence, we need only specify the external 
momentum relabelings.) Now numerator N 8 y should respect the 1 <H- 3 antisymmetry of 
the graph. This constraint, together with eq. (3.15), implies that 



N, 



S3' 



9 f 
S[U 

2 v 



t). 



The original cubic graph in fig. 11 then has the regulated numerator factor, 

9 



83 



■-k 4 s(u — t) 



(3.16) 



(3.17) 
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for the snail graph 83. Here k\ should not be set to zero until after this factor has canceled 
the 1/&4 propagator of the graph. The remaining two graphs, 84 and 85, can be determined 
in an analogous fashion using a nonplanar cut which also isolates a three-point subamplitude. 
Alternatively, the duality relations in appendix A fix them to be N 85 = N 8i = N 83 . 

We have carried out a very similar analysis of the corresponding cuts of M = 8 super- 
gravity and find that there are no snail contributions. This is in line with our heuristic 
expectations, described in section II D, that the gravity case should have an extra vanishing 
factor of kf in the numerator, setting all snail contributions to zero. 

Now that we have a complete Ansatz for the amplitude, the remaining task is to confirm 
that it has all desired properties and that it is a correct representation of the amplitude. 
Indeed, we verified that all numerators respect the graph symmetries and that all duality 
constraints hold on the sYM numerators. 

To prove that our construction is correct, we verified that a spanning set of .D- dimensional 
generalized unitarity cuts are properly reproduced. It has been shown [11] that there are no 
contributions to the four-loop four-point scattering amplitude in M = 4 sYM theory beyond 
those exposed on next-to-next-to-maximal cuts. This result was later confirmed using six- 
dimensional cuts [34]. As long as a candidate representation manifestly contains no worse 
loop-momentum behavior for individual terms than the previous representation, as is true 
of our current form, the set of next-to-next-to- maximal cuts are spanning (complete). We 
employ those cuts to verify our new expression, by comparing it against the previous one. 

We do not give the details here as the procedure is the same as given in refs. [11, 41], 
except that we generate the reference .D- dimensional analytic cuts using the previously- 
obtained forms of the amplitude [11]. For the gravity amplitude, obtained from the double- 
copy formula (2.5), we confirmed that its cuts through (next-to) 4 -maximal cuts match the 
corresponding cuts of the result in ref. [10]. 3 The agreement with the cuts of the earlier 
representation directly proves the duality and double-copy properties for the four-loop four- 
point amplitudes of N = 4 sYM theory and Af = 8 supergravity. 

IV. UV BEHAVIOR OF M = 8 SUPERGRAVITY AND M = 4 SYM THEORY 

In this section, we examine the UV properties of the four-loop four-particle M = 4 sYM 
theory and M = 8 supergravity amplitudes derived in the previous section, after reviewing 
lower-loop examples. Unlike the original form obtained for the M = 8 supergravity four-loop 
amplitude [10], the representation derived in this paper through eq. (2.5) is manifestly finite 
for D < 11/2. This property makes it much easier to determine its UV behavior. We will 
see that, in complete harmony with the corresponding M = 4 sYM amplitude, a divergence 
is indeed present in the expected critical dimension D c = 11/2. In M = 4 sYM theory this 
divergence in the amplitude corresponds to a counterterm of the schematic form Tr(T> 2 F 4 ), 
where T> represents a covariant derivative and F is the gauge field strength. Similarly, in 
M = 8 supergravity the amplitude divergence corresponds to a counterterm of the form 
V 8 R 4 , with the indices of the four Riemann tensors arranged in the supersymmetric combi- 
nation corresponding to the square of the Bel- Robinson tensor [65]. (For recent discussions 
of V 2k R A invariants in M = 8 supergravity see ref. [66].) We shall find a close connection 



3 The need for checking this high an order arises because of the large numbers of loop momenta occurring 
in individual terms in the original representation of the AT = 8 supergravity amplitude [10]. 
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between coefficients of the UV counterterms in these theories. 

In most of this section we will work in space-time dimension D = D c — 2e = 11/2 — 2e. 
In the last subsection we will discuss the UV behavior of the color double-trace terms in 
the M = 4 sYM amplitude at four loops, for which the divergence in D = 11/2 — 2e is not 
present. Here we move to the next dimension in which there can be a potential divergence, 
D = 6 — 2e, in order to assess whether there is one or not. 

In D = 11/2 — 2e, there are no UV divergences for either M = 4 sYM or M = 8 supergrav- 
ity below four loops; in both theories, the first divergences for L = 1,2,3 are for dimensions 
D c = 8, 7, 6, respectively. Therefore there can be no UV sub divergences in D = 11/2 — 2e at 
four loops. Near this dimension, the overall logarithmic UV divergences arise from integra- 
tion regions in which the loop momenta are much larger than the external momenta. For 
the representations of the M = 4 sYM and M = 8 supergravity amplitudes in this paper, 
each integral is either finite or develops at most a logarithmic divergence. We can capture 
the leading UV behavior by expanding in small external momenta [67, 68]. The logarithmic 
behavior of each integral means that we need only the leading term in this expansion. The 
amplitudes in question reduce to a collection of vacuum integrals that are relatively easy to 
evaluate. 

We begin by reviewing the UV properties of M = 4 sYM theory, in the light of the 
new representation of the four-loop four-point amplitude. Then we will turn to M = 8 
supergravity. 

A. Review of UV behavior of M = 4 sYM theory 

The UV properties of Af = 4 sYM theory at one through four loops were discussed in 
detail in ref. [11]. The first divergence appears in the critical dimension 

D c = 4 + i (£ = 2,3,4), (4.1) 

and D c — 8 at L — 1. In each case, the potential divergence is known to appear with a 
nonzero coefficient [11]. Thus the bound (4.1), proposed in refs. [8, 63], is in fact saturated, 
at least through four loops. 

For gauge group SU(N C ), gluon scattering amplitudes can be expressed in the trace basis, 
i.e. in terms of traces of products of generators in the fundamental representation. For four 
external gluons, only single-trace structures, of the form Tr (T ai T aj T ak T ai ), and double-trace 
structures, of the form Tr (T ai T a i) Tr(T ak T ai ) , can appear. At one and two loops, the color 
double-trace terms have the same UV behavior as the overall amplitudes, obeying eq. (4.1). 
However, starting at three loops they are less divergent [11]. (See also the discussions in 
refs. [64, 69].) In particular, the critical dimensions for finiteness for the double-trace terms 
satisfy 

D 2-trace = 4 § = 3, 4) . (4.2) 

Lj 

Using the representations of the three- and four-loop amplitudes described in previous sec- 
tions, we will see that these bounds follow more transparently than with the older repre- 
sentations in ref. [11]. That is because the color factors associated with the most divergent 
integrals now have a much simpler structure. 

Ref. [11] showed that the bound (4.2) is saturated at three loops. This computation did 
not involve any UV sub divergences because 4 + 8/3 = 20/3 < 7, and the first two-loop 
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y( Np ) 



FIG. 15: The two-loop vacuum integrals and 1/( NP ). Each (blue) dot on a propagator indicates 
an additional power of the propagator. 



divergences are at D c = 7. In section IV E we will show that the double-trace bound is 
saturated in the four-loop amplitude as well. In this case there are sub divergences, because 
4 + 8/4 = 6, and the first three-loop divergences are at D c = 6. 

Before proceeding to three and four loops, we review the two-loop case. The two-loop 
four-gluon amplitude as given in ref. [8] already obeys the color-kinematic duality. The full 
amplitude was originally presented in the trace basis [53]; in the critical dimension D c = 7, 
its divergence in terms of vacuum integrals is [11] 



A {2) 



SU{N C ) 
polo 



g 6 IC 



N 2 C V^ + 12(V^ + V^)') (4.3) 
x ( s ( Tri324 + Tri 423 ) + t ( Tri 24 3 + Tri 342 ) + u ( Tri 234 + Tri 432 



12 N c (V^ + V^ NP >) (s Tr 12 Tr 34 + t Tr 14 Tr 23 + u Tr 13 Tr 24 ) 



where and y( NP ), shown in fig. 15, are the fcf — > limit of the planar and nonplanar 
double-box integrals. 4 Each blue dot denotes an extra power of the propagator on which it 
lies. Here each dot coincides with the location of an external leg attachment in the original 
four-point integral; adjacent propagators separated by an external momentum become equal 
as that momentum vanishes. The factor 



K = sL4 tree (l,2,3,4) 



(4.4) 



contains the dependence on the external states. We use a shorthand notation for the color 
traces, 

Tr ijkl = Ti(T a *T a iT ak T ai ), Tr tJ = Ti{T ai T a ^) = 5 a ^ . (4.5) 
The poles of the vacuum integrals and y( NP ) in D = 7 — 2e are [8, 11], 



V (P) 

y(NP) 



7T 



20 (47r) 7 e 
30 (47r) 7 e 



(4.6) 



4 Wc normalize our integrals as in ref. [11], so that at two loops there is a relative minus sign compared to 
the normalization in ref. [81. 
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The UV divergence in the critical dimension, D c = 7, of the corresponding supergravity 
amplitude is [8] 



M 4 



(2) 



pole 



-2 (|) 6 stu(s 2 + t 2 + u 2 ) Mf ee (V^ + V^) . (4.7) 



We note that the UV divergence of the supergravity amplitude and that of the 1/N 2 - 
suppressed single-trace sYM amplitude are given by the same linear combination of vacuum 
integrals, namely + V^ NF \ We shall see that this pattern repeats itself through four 
loops. (This observation holds for the one-loop four-point amplitude as well. However, 
there the relation is rather trivial because both amplitudes are expressed in terms of the 
same scalar box integral; thus their UV divergences must be expressed in terms of a unique 
vacuum integral.) 

We now turn to three loops. The integrals that appear in the three-loop four-point M = 4 
sYM amplitude in the duality-satisfying form [30] all have 10 propagators. However, three of 
them, and I® in fig. 2, have one propagator depending solely on external momenta; 

we will refer to these integrals as nine-propagator integrals. From eqs. (2.15) and (2.7), their 
numerator factors are independent of the loop momenta: = AW = = (t — u)s/3. 

If the other integrals had two powers of the loop momentum in the numerator, then 
they would have the same generic large loop-momentum behavior as the nine-propagator 
integrals. However, the other integrals all have numerators that are at most linear in the 
loop momenta. Equation (2.19) shows that the numerator for the master graph (e) is 
linear in the r^, which are in turn linear in the loop momenta. The Jacobi relations (2.7) 
preserve this linearity for all other numerators. (In some cases the linear dependence cancels 
down to a constant behavior.) Therefore the leading UV divergence of the three- loop M = 4 
sYM amplitude comes from the three nine-propagator integrals. 

This result is consistent with a rearrangement of the leading UV terms of the earlier 
representation, as discussed in ref. [11]. Because graphs (j), (k) and (1) contain an external 
three-point tree, their color factors must be proportional to the product of two structure 
constants, multiplied by a color Casimir operator which has no free color indices [11]. The 
product of the two structure constants, j a ^ h j ha '-i a * or a permutation thereof, takes the same 
form as a tree amplitude, namely a single color trace. Hence the leading UV divergence at 
three loops, in the critical dimension D c = 4 + 6/3 = 6, contains no double color-trace 
terms [11], while they are present in the leading divergence at two loops, eq. (4.6). The 
duality-satisfying representation of the amplitude automatically has the no-leading-double- 
trace feature. Very similar behavior is observed within a string theory analysis [64]. There 
the leading behavior at three (and four) loops is dominated by the collision of pairs of vertex 
operators, producing inverse momentum factors, reminiscent of the form of graphs (j), (k) 
and (1). The finiteness of the three- loop double-trace terms in D = 6 remains puzzling from 
the point of view of field-theoretic algebraic nonrenormalization considerations [69]. 

To extract the UV divergence, we carry out the small momentum expansion [67, 68]. 
Integral (j) reduces to the vacuum integral [11, 55] displayed in fig. 16, while integrals 
(k) and (1) reduce to . These integrals diverge first in D = 6. Their color factors 
are closely related: We can use a color Jacobi identity involving the boxes in the upper 
right-hand corner of graphs (j) and (k) of fig. 2 to show that the difference between Cq) and 
C(k) contains a triangle subgraph. A further color Jacobi identity allows us to replace the 
triangle by a three-vertex, multiplied by the quadratic Casimir factor Ca/2 = N c . Iterating 
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FIG. 16: The three-loop vacuum integrals and TA B ). 



this procedure, and also applying it to the difference of Cm and Cm, we find that 

Cq) - N 3 J a ^ b f ba ^ = C (k) = C (1) . (4.8) 

Note that C(k) and Cm are associated with nonplanar graphs, and hence have only 
subleading-color terms. Equation (4.8) states that the subleading-color parts of Cm, Cha 
and C(i) are all equal. 

Taking into account that the combinatorial factor of I® is twice as large as those for 7^ 
and 7^ k ^ in fig. 2, and expressing the color factors for (j), (k) and (1) in the trace basis, it is 
straightforward to see that the UV divergence in D c = 6 is [11], 
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The UV poles of the vacuum graphs are [55] 
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(4.11) 



We may compare this result to the UV divergence in D c = 6 of the corresponding M = 8 
supergravity amplitude 



pole 



[10(^ A ) + 3V^)] 



(4.12) 



In this case the leading UV divergence involves not only the 1PR nine-propagator integrals, 
but also the 1PI ten-propagator integrals whose sYM numerators are linear in the loop 
momenta, because their (squared) supergravity numerators are quadratic. Once again, as 
was true for the two-loop amplitude, the UV divergence of the three-loop supergravity 
amplitude and that of the 1/iV^-suppressed single-color-trace three- loop sYM amplitude are 
given by the same combination of vacuum integrals, namely + 3 . Below we shall 
see that the same phenomenon persists through four loops. 

Using the representation of the four-loop four-point M = 4 sYM amplitude derived in 
earlier sections it is equally straightforward to recover its UV divergence. Indeed, an inspec- 
tion of the integrals listed in figs. 5-12 and of their numerator factors listed in appendix B 
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FIG. 17: The basic four- loop vacuum integrals V±, V2 an d to which all others can be reduced. 

reveals that the leading UV behavior comes solely from integrals /go through 7 85 . These six 
integrals have 11 internal propagators, and numerator factors that are independent of the 
loop momentum. Therefore they diverge first in D c = 11/2, which matches the expected 
critical dimension, D c = 4 + 6/L with L = 4. 

In contrast, the 1PI integrals I\ through J52 have 13 internal propagators. Their nu- 
merators would have to be quartic in the loop momenta for them to diverge in D = 11/2. 
However, we note from eq. (3.14) that the master numerators N±s and A 2 g are quadratic in 
the Tij, and hence merely quadratic in the loop momenta. The Jacobi relations (A4) and 
(A5) preserve this quadratic behavior for all numerators. Therefore integrals l\ through J 52 
are finite in D c = 11/2. The 1PR integrals J53 through I? 9 have 12 internal propagators. If 
their numerators were quadratic in the loop momenta, then they would diverge in D = 11/2. 
However, it is easy to see from eqs. (A5) and (Bl) that their numerators are all linear in the 
loop momenta. 

Integrals /go through Jgs reduce easily to vacuum integrals in the limit that the external 
momenta vanish. The planar integrals /go and Jg3 reduce to the vacuum integral V\ depicted 
in fig. 17. While integrals igi and I^a are nonplanar as four-point graphs, in the vacuum 
limit they reduce to the planar vacuum integral V2. Finally, integrals /§2 and reduce to 
the nonplanar vacuum integral Vg. 

As was the case at three loops, the color factors for the leading UV graphs are related 
by color Jacobi identities. In this case we can subtract, for example, Cgi from Cgo, and use 
a Jacobi identity operating on the box at the top center of the graphs in fig. 11. Then we 
reduce the resulting triangle subgraphs iteratively, to find 



Again, the subleading-color parts of all contributing color factors are equal. 

From eq. (A5), the numerator factors obey A 80 = Agi = A 82 and A 83 = A r 84 = A^ 85 . 
On the other hand, the combinatorial factor of igi in eq. (3.3) is twice as large as those 
for J 80 and / 8 2, and similarly for / 84 with respect to /g3 and / 8 5. Taking into account both 
combinatorial and numerator factors, the contribution of / 83 is — ~ times that of /go, and 
similarly for the other two pairs of graphs. Combining all terms and switching to the color- 
trace basis, we find that the UV divergence in the critical dimension D = 11/2 is given 
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in agreement with the results of ref. [11]. It is interesting to note from eq. (4.14) that the 
single-trace UV divergence in D = 11/2 has and N% components, but the N® component 
vanishes. There is currently no general explanation for this fact, apart from the explicit 
values of the color factors C 8 o through C 85 . 

The values of the three master integrals appearing in eq. (4.14) are [11] 
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where NO m denotes a certain three-loop two-point nonplanar integral. While its analytic 
expression in D = 11/2 is not known, it may be evaluated numerically [11] using the Gegen- 
bauer polynomial x-space technique (GPXT) [70], with the result NO m = —6.1983992267 . . .. 
We remark that an evaluation to much higher accuracy is also possible [71] using the DRA 
method [72] , which involves combining dimensional recurrence relations [73] with analyticity 
in the space-time dimension. (This latter method has been applied to similar, and even more 
complex, integrals in refs. [74].) 

In the following two subsections we will obtain the UV divergence of the corresponding 
M = 8 supergravity amplitude. We will find that, as for the two- and three-loop case, it is 
given by the same combination of vacuum integrals as the 1/iV^-suppressed single-trace UV 
divergence of the corresponding four-point M = 4 sYM amplitude, namely V\ + 2 V% + V%. 



B. M = 8 supergravity vacuum graphs at four loops 



With the duality-satisfying form of the corresponding M = 4 sYM amplitude as the start- 
ing point, the double-copy formula (2.5) immediately gives us an expression for the four-loop 
M = 8 supergravity integrand. We confirmed this integrand by comparing its cuts with the 
cuts of the known four-loop Af = 8 amplitude [10] . The double-copy formula is equivalent 
to the following squaring relation for the M = 8 supergravity numerators: 

N M=8 = (ivf=4)2_ ( 416 ) 

By inspecting the squares of the numerator factors listed in appendix B, and counting the 
number of loop momenta in the numerator of each integral, it is easy to see that, in all cases, 
the integrals composing the resulting M = 8 supergravity amplitude are manifestly finite for 
D < 11/2. 

Indeed, as remarked in the previous subsection, the maximum degree in loop momenta 
of the numerator factors of the sYM amplitude is 2 for the 13-propagator integrals, 1 for 
the 12-propagator integrals and for the 11-propagator integrals, where we count only 
those propagators carrying loop momentum. Consequently, the maximum degree in the 
loop momenta of each supergravity numerator polynomial is 4, 2 and for the 13-, 12- and 
11-propagator integrals, respectively. Such integrals all generically diverge logarithmically 
in D = 11/2. Thus, the worst UV behavior of any M = 8 supergravity integral matches that 
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of the worst-behaved integrals for the J\f = 4 sYM amplitude, i.e. the 11-propagator graphs 
(80)-(85) of fig. 11, which have the form of propagator corrections. The main difference 
is that now graphs other than propagator corrections carry this leading behavior. The 
representation of the four-loop amplitude described here reproduces the finiteness bounds 
of ref. [10], but the UV behavior is now manifest, allowing us to avoid performing any loop 
integration to expose this feature. 

The new representation also makes it far simpler to determine whether the finiteness 
bound D < 11/2 is saturated, by extracting the precise UV divergence of M = 8 supergravity 
in D c — 11/2. In our earlier representation [10], the UV divergence required the sixth-order 
terms in the expansion in small external momenta, making it rather cumbersome to extract. 
Now that the UV behavior is manifest, only the leading term in the expansion is required. 
This feature means that for each integral, we need only retain the terms in each sYM 
numerator with the highest powers of the loop momenta, and then square them. 

The result of the expansion in the external momenta is a collection of tensor integrals, 
in which the numerator factors have a homogeneous degree in the loop momenta and are 
polynomials in the scalar products of loop and external momenta. Such integrals may be 
further reduced by making use of Lorentz invariance in order to extract the dependence on 
the external momenta from the tensor integrals. More precisely, under integration we can 
replace a generic two-tensor by 

H- ^V^k-lj, (4.17) 

and a four-tensor by 

l^ifM^m ^ (Ar)^ r)^ 1 +Br]^ k 77 ww + Ct? ww f») , (4.18) 

J [D — 1)D[D + 2) 

where 

A = (D + l)/j ■ lj Ik ■ k — h - h lj ■ k — h • k lj ' h , 

B = —li ■ lj Ik ■ h + (D + ■ Ik lj ■ k — h ■ h lj ■ h ? 

C = -k ■ lj Ik-k-k- h lj -k + (D + l)k ■ k lj ■ l k . (4.19) 

Upon using these identities and summing over all permutations we find that the UV pole 
has the form, 



M 4 



_I (1) W s tu (s 2 + t 2 + u 2 ) 2 Mf ee V (4) , (4.20) 



pole 4 V 2 

where M^ ee is the four-point tree-level supergravity amplitude, and 
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V (4) = ^ (4.21) 



is the sum of integrals 7" shown in figs. 18-20. All the kinematic dependence has been 
extracted in eq. (4.20); each integral IV is a pure number multiplied by a 1/e pole. The 
numerator factors for these integrals are given in the first column of table I (not counting the 
column labeling the integrals IV) in appendix C. The overall kinematic factor of (s 2 +t 2 +u 2 ) 2 
is guaranteed by complete permutation symmetry of the four-point amplitude. 
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(1) (2) (3) (4) (5) 




(6) (7) (8) (9) (10) 




(11) (12) (13) (14) (15) 




(16) (17) (18) (19) (20) 




(21) (22) (23) (24) (25) 

FIG. 18: Vacuum graphs If through i^s- The momentum labels refer to the internal lines carrying 
an arrow. Their numerator factors are listed in table I in appendix C. 

The precise combination of vacuum integrals 1%, in particular their numerical coefficients 
and numerator tensor structures, is sensitive to the choice of independent loop momenta 
in each of the integrals appearing in the complete supergravity amplitude. This choice is 
inherited in turn from the parametrization of the integrals of the M = 4 sYM amplitude, 
through the squaring relation (4.16). The ambiguity in choosing the loop momenta may also 
be used to generate identities between different integrals. 

We now discuss in some more detail how we arrived at the results in the first column 
of table I in appendix C. After applying eqs. (4.17) and (4.18), the numerator factors of 
the vacuum integrals appearing in which no longer carry any external momenta, can 
be simplified further using loop-momentum conservation. Often, inverse propagators can be 
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FIG. 19: Vacuum graphs I| 6 through I§ Q . 



identified, using for example U-lj = If— I]) when k+lj = Ik corresponds to a propagator 
for the graph. Such factors will cancel existing propagators and lead to simpler integrals 
with fewer propagators and lower-degree numerator factors. In carrying out this procedure, 
it is useful to ensure that propagators are not 'over-canceled'. For example, no squared 
inverse propagator should appear in a vacuum integral that does not have the corresponding 
propagator raised to at least the second power. For certain integrals the complete numerator 
factor may be written as a combination of different inverse propagators. In this case, the 
integral reduces to a combination of scalar integrals of different topologies, each of which is 
obtained from the initial integral by collapsing some of its propagators. For other integrals, 
the rank of the original numerator tensor is reduced by only two units. For 13 integrals the 
numerator factor remains of fourth order in loop momenta with no canceled propagators; in 
all such cases the numerator factors carrying loop momenta may be expressed as a perfect 
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square, (ij • Z,) 2 for different z and j. 

To illustrate this procedure in more detail, we consider two examples — the reduction of 
integrals I§§ and As we will see, the former integral reduces to only scalar integrals, while 
the latter integral leaves behind a four-tensor integral. The numerator factor of the integral 
Iqq in the supergravity amplitude is given by the squaring relation (4.16), and formula (Bl) 



for A^= 4 , to be 



N: 



JV=8 



4t(r 35 - 2r 36 ) + 2m(t 35 + 3r 45 - 4r 46 ) - s(6u + r 1B - 6* + 5r 25 - 8r 26 ) 



(4.22) 



For this 12-propagator integral, the leading UV behavior comes from terms that are quadratic 
in the loop momenta in N^ =s , so the two r-independent terms inside the brackets in 
eq. (4.22) may be dropped. For the remaining terms, eq. (4.17) implies that iV 66 is equivalent 
under integration to 



-N; 



Af=8 



66 



H> — s 2 tu(llT 55 - 48r 56 + 48r 66 ) - — s 4 (15r 55 - 64r 56 



64r, 



66 I 



(4.23) 



We have extracted a factor of 1/s on the left-hand side from the external propagator in the 
1PR integral I 66 . 
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Next we express all scalar products of internal momenta in terms of inverse propagators, 
using r 55 = 2l\ and tqq — t$q = —l\ + 1\ + (I5 — Iq) 2 . Summing over the 4! permutations, and 
dividing by the symmetry factor of 4 for this graph, we find that the contribution of Iqq is 

;< =8 H- ^( S 2 + t 2 + n 2 ) 2 (17/ 5 2 - 32(/ 5 -/ 6 ) 2 - 32^). (4.24) 

In the permutation sum the first term in eq. (4.23) cancels out completely due to s+t+u = 0. 

Each one of the factors I 2 , (1 5 — Iq) 2 and l\ in eq. (4.24) cancels one propagator of 
The topology resulting from l\ can be identified as graph (17) in fig. 18. The (I5 — Iq) 2 and 
Iq terms both lead to graph (1) in fig. 18. Therefore, to leading order in the small external 
momentum expansion in D = 11/2, the contribution of integral Iqq to becomes 

J 66 ^ g(17ir 7 -64/r) = ^Ilr-^P,. (4.25) 

Here the I? denote the vacuum integrals with topologies shown in figures 18-20, and 
momentum-dependent numerator factors shown in the first column of table I, but with 
the rational numerical coefficients set to unity. For example, with this notation we have 

= ~ ^Iss^ i • ^ e ^ ave mm tiplied by 4/(s 2 + t 2 + w 2 ) 2 in passing from eq. (4.24) to 
eq. (4.25), to account for the relative prefactor of (s 2 + t 2 + u 2 ) 2 /A in eq. (4.20), compared 
with eq. (3.2). The difference between the numerical coefficients in eq. (4.25) and those in 
table I is simply that the table collects the contributions to each vacuum graph from all the 
different integrals li appearing in eq. (3.2). 

As a second example, consider the reduction of 1 12- The numerator factor in the M = 8 
supergravity amplitude is 



M=8 b 



s (ri6 - r 26 - t 35 + t 45 + 2t) 

1 2 

+ 2 s r 56 - 2 (4ti 6 t 25 + 4ti 5 t 26 + r 45 (r 36 - 3r 46 ) + T 35 (r m - 3r 36 )) . (4.26) 

The leading UV divergence requires four powers of the loop momenta from the numerator. 
That can come only from the second and third terms inside the square; the first term is linear 
in the loop momentum and can be dropped. Using eq. (4.18) and the on-shell condition for 
the external legs we find, after some straightforward calculation, that under the integral sign 
A^i2 =8 is equivalent to 

^ ~ (D - 1)D(D + 2) sHu " 2)T ™ + jDr55T66 ] (4 ' 27) 

+ {D-l)D{D + 2) sA [ (j ° 3 ~ 19j ° 2 + U6D ~ 96)T ™ + 4(17j ° " 25 ) T55T66 ] ' 

By inspecting the graph for I12 it is straightforward to see that cannot be completely 
expressed in terms of inverse propagators. We therefore keep it in this form. The sum over 
the 4! permutations leads, as in the case of A<g^ =8 , to the complete cancellation of the term 
proportional to s 2 tu. Dividing by the symmetry factor of 4 for 7 12 , we find that 



A^ =8 ' — ^ (s 2 + t 2 + u 2 ) 2 



D 6 — 19-D 2 + 146-D — 96 o V7D - 25 ,„,, 
! T- 2 -1-16 / / 

(D - 1)D(D + 2) 56 (D — 1)D(D + 2) 5 6 



. (4.28) 
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Considering the graph associated with I12, we see that the r| 6 term generates graph (3) in 
fig. 18, whereas the factor of /|/g reduces the tripled propagators to doubled ones, generating 
graph (4) in fig. 18. Multiplying by an overall factor of 4/(s 2 +t 2 +u 2 ) 2 and setting D = 11/2, 
we obtain for the contribution of I12 to V*- 4 ^, 

9556 - 35072 - 
/l2 H> 1485 /3 + U85 h ' (429) 

It turns out that integral ii 2 generates the only contribution to J3 , so that Jg = |||| ^3 , as 

given in table I in appendix C. Several other integrals li contribute to 1%, so the rational 
coefficient given in table I for l\ is different from the one in eq. (4.29). 

Carrying out the steps detailed above for all the integrals appearing in the amplitude, we 
find that they reduce to 92 different vacuum integral topologies. In the complete amplitude 
some of them cancel out 5 , leaving only the 69 vacuum topologies shown in figures 18-20. 
While these integrals may be rearranged somewhat using momentum conservation, to find 
the full set of relations between integrals we need more powerful methods. We shall see in 
the next section that integration-by-parts identities reduce them all to linear combinations 
of the three scalar vacuum integrals shown in fig. 17. 



C. Integrating the vacuum graphs 

To evaluate vacuum integrals in their critical dimension D c = 11/2, we use the same 
infrared rearrangement [67] (related to the R* operation [75]) that we used previously to 
evaluate the three-loop vacuum integrals and [55] and the four-loop vacuum inte- 
grals Vi, V2 and Vg [11]. This approach was already discussed in some detail in ref. [11], so 
here we include only a brief summary. 

The expansion in small external momentum introduces unphysical infrared divergences. 
To separate them from the UV singularities we inject and remove momentum k^, with k 2 ^ 0, 
at two of the vertices of the vacuum integral, thus transforming it into a four-loop two-point 
integral. This two-point integral possesses the same UV poles as the vacuum integral, but no 
infrared divergences. We always take the two vertices in question to be connected by a single 
propagator in the four-loop vacuum integral. In this case, the four-loop two-point integrals 
can be factorized into products of three-loop two-point integrals and one-loop two-point 
integrals. The one-loop two-point integrals are trivial, and contain the UV 1/e pole. The 
finite three-loop two-point propagator integrals are then evaluated through the method of 
integration by parts (IBP) [76], implemented in the MINCER algorithm. For certain tensor 
integrals we also applied the Laporta algorithm [60] for solving the integration-by-parts 
relations, as implemented in the computer code AIR [77]. In addition, we employed gluing 
relations [76], which demand consistency of the various ways of factorizing the four-loop 
UV-divergent integral into products of lower-loop integrals. These consistency conditions 
are nontrivial and aid in the evaluation of some of the master integrals remaining after IBP 
reduction in D = 11/2. 

Further consistency relations between vacuum integrals can be derived as follows. We 
start with a four-loop two-point integral with a numerator factor of degree less than or 
equal to two in the loop momenta. We expand it to next-to-next-to-leading order in the 



5 Some of these cancellations are dimension-independent, while others occur only in D = 11/2. 
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small external momenta, using different parametrizations of the loop momenta. We require 
that the different parametrizations yield consistent results. These relations may be used 
either to reduce the number of integrals that need to be evaluated, or as consistency checks 
of the integral evaluation. A similar strategy was used to show that the four-loop four- 
point M = 8 supergravity amplitude is UV finite in D = 5 [10]. In that case, the complete 
amplitude was expanded for different loop-momentum parametrizations. This approach led 
to consistency relations involving only integrals that naturally appear in the amplitude. The 
present situation is somewhat more involved because evaluation of each of the 69 integrals 
implies that we need to expand integrals which do not necessarily appear in the final result. 
In this case, the consistency relations that are generated involve a much larger set of integrals 
than just the 69 under consideration. 

Following the strategies reviewed above, all the 69 integrals appearing in the small mo- 
mentum expansion of the four-loop four-point M = 8 supergravity amplitude can be reduced 
to linear combinations of the integrals V\, V2 and V§ shown in fig. 17. These integrals have 
already appeared in the UV divergence of the corresponding M = 4 sYM amplitude. The 
coefficients of Vi, V 2 and V§ after this reduction are provided in the second, third and fourth 
columns, respectively, of table I in appendix C. The values of the UV poles of V\, V2 and 
V$ were determined in ref. [11] and are given in eq. (4.15). 

It is interesting to complete the evaluation of the two examples discussed in the previous 
subsection, 7 66 and J12, and thereby illustrate a general feature regarding the positivity prop- 
erties of the residues of the UV poles of the integrals appearing in the Af = 8 supergravity 
amplitude in the double-copy representation. 

The results of appendix C can be used (after dividing by rational coefficients in the second 
column) to express the following integrals 7" in terms of the basis integrals V±, V2 and Vg: 

Ti = V 2 , (4.30) 

~ = 6671 _ 531 

3 4800 1 800 2 ' K ' 1 

It = V l , (4.32) 

^7 = ^1 + ^2, (4.33) 

Inserting these equations into eq. (4.25), for the contribution of I§§ to V < - 4 - ) we obtain 

Ja6 ^ 544 ^ 2048 - = 408 1776 

66 11 17 11 1 11 11 V 1 

Using the expressions (4.15) for V\ and V 2 we find that the sign of the UV divergence of 7 66 
is negative. 

Similarly, the contribution of the UV singularity of I12 to is 

9556 f 35072 f 58023419 MOggl 

1485 3 1485 4 1782000 33000 v ' 

Inserting the numerical values of V\ and V 2 into eq. (4.35) we find that, unlike the UV 
pole of 7i2 has a positive residue in D = 11/2. 

We see that the residue of the leading UV pole can be either positive or negative, despite 
the fact that the initial numerator factors are perfect squares. This phenomenon occurs 
frequently in the reduction of the various integrals appearing in the N = 8 supergravity 
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amplitudes. It leads to strong cancellations among the corresponding UV poles. The origin 
of this phenomenon is simply that the propagators occurring in the denominators of the 
integrand do not have a fixed sign. For example, the 1PR integral Iqq contains an explicit 
factor of 1/s from a propagator external to the loops. The sum over its permutations contains 
factors of 1/t and 1/u. Not all of these factors can be positive, given that s + t + u = 0. 
Even if these factors are multiplied by positive numbers, the sum can be negative. Indeed, 
in any physical region (with one Mandelstam invariant positive and two negative), if they 
were all multiplied by the same number, the sum would be negative, because - + 7 + - < 
for s + t + u = 0. 



D. The M = 8 supergravity UV behavior 

The expressions for the leading UV divergences of the 69 vacuum integrals IV in D = 
11/2— 2e are collected in table I in appendix C. Totaling up these contributions, we determine 
the value of defined through eqs. (4.20) and (4.21), to be 

~ 93 

V (4) = -f (Yi + 2V 2 + V 8 ) . (4.36) 

Plugging this value for into eq. (4.20), the UV singularity of the four- loop four-point 
Af = 8 supergravity amplitude is given by 



M 



(i) 



23 /K\i° 

pole 8 V 2 



(_j s tu(s 2 + t 2 + u 2 ) 2 Ml Tee (V 1 + 2V 2 + V 8 ). (4.37) 



Because this result is nonvanishing, the four-loop four-point amplitude of Af = 8 supergravity 
diverges in exactly the same critical dimension D c = 11/2 as that for Af = 4 sYM theory. 

Comparing eqs. (4.37) and (4.14), we see that the leading UV divergence of the four- 
loop four-point Af = 8 supergravity amplitude is given by the same combination of vacuum 
integrals as the l/A^-suppressed single-trace UV divergence of the four-point Af = 4 sYM 
amplitude, namely V\ + 2V 2 + V s . This pattern matches our observations above regarding 
the two- and three-loop amplitudes. While we do not understand the full significance of 
this close relation between the UV divergences of the gravity and subleading-color sYM 
amplitudes, it seems unlikely to be accidental. If it persists at higher loops as well, it 
could potentially have interesting consequences for the higher-loop UV behavior of Af = 8 
supergravity, especially given that Af = 4 sYM is UV finite in four dimensions. The fact 
that this nontrivial connection holds through at least four loops is a central observation of 
this paper. 

While the divergences of all integrals entering A4^ are expressed in terms of the vacuum 
master integrals V±, V 2 and Vg, the expression (4.37) arises as result of extensive cancellations 
between the various integrals. For example, integrals /go through 7 85 are responsible for the 
complete UV divergence of the four-loop four-point A" = 4 sYM amplitude. It is easy to 
see that they contribute to the UV divergence of the Af = 8 supergravity amplitude the 
following amount: 



Mi 4) 



64 (-J stu (s 2 + t 2 + u 2 ) 2 cc (Fx + 2V 2 + V s ) . (4.38) 



^80,81,82 _ _ /K\ 10 

pole 



Note that the snail integrals Jg3, 7g4 and I$5 yield vanishing contributions in the supergravity 
case. Because the rational coefficient in eq. (4.37), 23/8, is smaller than 64, it follows that 
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FIG. 21: The vacuum integrals U±, U2 and U3 that would be generated by the graphs 51, 52 and 72 
in the double-trace part of the four-loop divergence of W = 4 sYM theory in D = 6. Because these 
three graphs possess both an overall divergence and a three-loop subdivergence, we only evaluate 
them after first performing a subtraction of the subdivergence. 



the integrals I\ through I 79 contribute to the UV pole with opposite sign, relative to I 80 , hi 
and I$2- This fact is not surprising: although the numerators are squares, the propagators 
do not have all the same signs. 

To expose a more dramatic numerical cancellation, we have split the above computation 
into two pieces: The 12-propagator contributions, denned to be those from integrals J53 
through J79, except J72, which all have explicit single inverse powers of the Mandelstam 
variables, s,t or u; and the remaining 11- and 13-propagator contributions, which have 
either no such explicit factor, or else the factor is squared. We find that 



-Mi 4) 



12-propagator 



pole 

1 1 , 13-propagator 
pole 



- 142 (2 

1159 



10 



stu (V 



t z + u 2 ) 2 Mf ee (Fx 



2V 2 + V 8 



10 



stu (s 2 + 1 2 



2\2 



U 



Mf ee (vi + 2v 2 + v 8 . 



(4.39) 
(4.40) 



Remarkably, the numerical cancellation between these two sets of contributions, in order to 
get the total (4.37), is quite significant, to within about 2%. That is, (23/8) /142 = 0.0202 . . .. 

In summary, by carrying out the integration to extract the four-loop four-point UV 
divergence of M = 8 supergravity in D = 11/2, we have found extensive cancellations, both 
analytical and numerical, once again revealing surprising hidden structure. 



E. The M = 4 sYM color double-trace UV divergence 

In ref. [11] we showed explicitly that the color double-trace terms in the four-point M = 4 
sYM amplitudes are better-behaved in the UV than the single-trace terms, starting at three 
loops. In particular, they obey the finiteness bound D < 4 + 8/L ah three and four loops. 
This improved behavior has been discussed from the vantage points of both string theory 
and field theory [64, 69]. The form of the four-loop four-point sYM amplitude constructed 
here allows us to evaluate the four-loop double-trace divergence in D = 6 — 2e dimensions 
and probe whether the bound is saturated. We will find that the double-trace term does 
indeed diverge, and so the bound is saturated. A subtlety arises in this calculation because 
the single-trace terms in the three-loop four-point amplitude diverge in six dimensions. 
Therefore the extraction of the four-loop four-point double-trace divergence requires a careful 
subtraction of the contribution of the three-loop counterterm. 

In light of the presence of sub divergences in D = 6 — 2e, it is useful to separate the 
integrals appearing in the four-loop amplitude into those that cannot have sub divergences 
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FIG. 22: Three-loop propagator subgraphs P5i(Z|), -P52G5. 
in graphs 51, 52 and 72, respectively. 




P^Jlh 

and -P72G5), generated by setting fc 4 — > 



(graphs 1 through 50) and those that might (graphs 51 through 85). The extraction of UV 
divergences uses precisely the same methods as discussed above in the analysis of the UV 
behavior of the M = 8 supergravity amplitude. For the first set of graphs, we simply quote 
the result without providing further details. The UV divergence from graphs 1 through 50 
has the following Tri2 Tr34 component: 



4 4) 



Tr 12 Tr 34 
pole 1-50 



9 10 IC 



3(4vr; 
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Tr i2 Tr s4 N, 



31 



; s 2 + f 2 + u 2) ( N 2 {1 _ + 10Cs) + 180Cs ) _ 9,2 ^2 Ca + 25C ^ ^ (441) 



where K is defined in eq. (4.4). The other double-trace components can be obtained by 
permuting the external momenta in this expression. The absence of a 1 /e 2 pole signals that 
this class of graphs indeed has no subdivergences. 

Next we evaluate the remaining graphs, 51 through 85. Because the color factors of 1PR 
graphs do not contain double traces, it follows that, in fact, only the integrals /51, -Z52 and I72 
contribute to the double-trace terms. If they did not have subdivergences, we could evaluate 
them directly by setting the external momenta to zero, leading to the vacuum integrals U\, 
U2 and U3 shown in fig. 21. Instead, we will perform a subtraction of the subdivergence and 
evaluate the inner three- loop integral first, before evaluating the integral over the outer loop 
momentum. 

We first note from eq. (A5) that the numerator factors of these three integrals are all equal, 
7V51 = N52 = N-J2- All three integrals contain an essentially identical subdivergence, from a 
three-point three- loop subgraph whose external legs carry momentum k^, l§ and I5 + k^. In 
the evaluation of this inner graph, and its subtraction term, if we are only interested in the 
final 1/e 2 and 1/e contributions, we can neglect the dependence on /c 4 : the 1/e 2 divergence is 
independent of the details of the external momenta, and the 1/e contribution, arising from 
the finite part of the three-loop subgraph, comes from the integration region Z 5 ^> k& where 
the outer loop diverges. Then the three-point subgraph reduces to a propagator (two-point) 
subgraph, as shown in fig. 22 for the respective cases of graphs 51, 52 and 72. The blue dot 
indicates the location of the doubled propagator that is generated in the limit k^ — > 0. 

These integrals can be evaluated in D = 6 — 2e using IBP identities and gluing relations 
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through the necessary order, 0(e°). The results are: 
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(4.42) 
(4.43) 
(4.44) 



We note that the propagator subgraphs for graphs 52 and 72 happen to have identical values 
through 0(e). Also, inspecting the form of the integrals, and comparing with eqs. (4.10) 
and (4.11), we see that the leading 1/e singularities correspond to the three-loop vacuum 
integrals (for the planar graph 51) and (for graphs 52 and 72). 

In the trace basis, the color factors for graphs 51, 52 and 72 have the following form: 



051 



0] 



52 



iV 2 (iV 2 + 12)(Tr 1234 + Tr 1432 ) 

+ 2iV c (iV 2 + 12) ( Tr 12 Tr 34 + Tr 13 Tr 24 + Tr 14 Tr 23 ) , (4.45) 
C 72 = 12iV c 2 (Tr 1234 + Tr 1432 ) + 24iV c ( Tr 12 Tr 34 + Tr 13 Tr 24 + Tr 14 Tr 23 ). (4.46) 



The double-trace parts of these color factors contain a piece proportional to iV 3 and one 
proportional to iV*. The iV 3 part comes only from graph 51. The parts of these color 
factors are equal for the three graphs. Taking into account the relative symmetry factors 
in eq. (3.3), we see that the relevant linear combination of propagator integrals for the 
part is P N i = P51 + 2P 52 + P72, which is given by, 



Pm{l 
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-3-ye 



(4tt; 
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f C 5 + 0(e) 



(4.47) 



Next we need to identify a subtraction that accounts for the three-loop counterterm 
needed to cancel the pole given in eq. (4.9). We could compute the four-loop counter- 
amplitude by sewing a tree amplitude onto the matrix element of the counterterm (essentially 
the negative of eq. (4.9)), in close analogy to how we evaluated the UV divergence of the 
three-loop four-point supergravity amplitude in odd dimensions above D = 6 [55] . However, 
it is safer to perform the subtraction within the integrals for the individual graphs 51, 52 
and 72, in order to ensure that there are no spurious contributions arising from potentially 
different ways of regularizing the infrared behavior. 

In discussing higher-loop divergences of a theory, there is always a freedom associated 
with additional, finite renormalizations of the theory at lower loops, in this case three loops. 
Here we will choose, for definiteness and simplicity, an MS scheme for the three-loop renor- 
malization. In this scheme, the necessary counterterms are 
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Notice that the factor of (— 1\) 
contributions (4.48) and (4.49). 



(4vr) 9 - 3e 2e ' 

" 3<E in eqs. (4.42) and (4.47) is absent in the counterterm 
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The evaluations of the remaining one-loop integrals and of their subtraction are essentially 
identical for the graphs 51, 52 and 72. We will discuss in detail only the graph 51 for which 
we consider the following subtracted integral, 
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(4.50) 



If a term in the outer integral over l 5 is UV convergent, one can set e — > and the 1/e 
poles will cancel in the brackets, leaving no UV contribution at all at four loops. On the 
other hand, if a term in the outer integral generates a logarithmic divergence, as in the 
terms in N 5 i that are quadratic in Z 5 , then the leading 1/e 2 pole from -Psi(Zf) is not cancelled 
by the counterterm contribution P§i'. Accounting for a factor of (— l 2 )~ e from the one-loop 
integration measure, we see that there is a mismatch by a factor of four. That is, the leading 
1/e 2 contribution from Pg£ m is of opposite sign to that from P 51 (Z 2 ) and is four times larger 
in magnitude, due to a factor of (— l 2 )~ e instead of (— /f) _4e . 

The relative factor of four originates physically from the locality of the UV poles in 
the subtracted four-loop amplitude. The four-loop amplitude carries fractional dimension 
oc s _4e , while the one- loop counterterm amplitude carries fractional dimension oc s~ e . Ex- 
panding these s-dependent factors in e, we see that a ratio of (—4) between the 1/e 2 poles 
from the amplitude and its counter-amplitude is required, in order for non-local terms of 
the form 1/e x ln(— Sy) to cancel in the subtracted amplitude. 

From eq. (Bl), the quadratic terms in the numerator factor A^i are given by, 



We insert eq. (4.51) into eq. (4.50) and Feynman parametrize the resulting one-loop integral. 
As usual with a one-loop integral with numerator quadratic in the loop momentum, there 
are two types of terms: 

1 . terms that depend in a fairly complicated way on the Feynman parameters and external 
momentum invariants, originating from the shift l§ = + A M in the loop momentum 
Z5 needed to complete the square in the denominator, and 

2. terms from integrating over the shifted loop momentum g M , in which one can use the 
identity q^q v = q 2 / D x rf v . 

The former terms do not contain an ultraviolet divergence from the outer integral; therefore 
they can be dropped as discussed above. In the latter terms, the first set of terms in 
N^ ad in eq. (4.51), containing the —6 prefactor, drop out because k\ = k\ = k\ = 0. In 
the second set of terms in A r g| iad , one can effectively make the replacement (4.17), because 
any extra terms due to external momentum dependence are finite and unrelated to the UV 
divergence and therefore drop out. (We checked this statement by performing a full Feynman 
parametrization.) 

Therefore, the 1/e 2 and 1/e terms in I^f 3 in eq. (4.50) are correctly captured by 
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We do not need to evaluate this full integral; to extract the UV pole it is sufficient to simplify 
it to the form of a massive bubble integral by rearranging the external momenta, 
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In the term coming from P§i(l 2 ) in eq. (4.53), the powers to which the two propagators in 
the D = 6 — 2e bubble (after analytic continuation) are raised are 2 + 3e and 1; while in the 
term coming from P^' they are 2 and 1. 

The momentum-independent parts of these bubble integrals are given by 
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(Other ways of regulating the infrared behavior, such as keeping the complete dependence 
on external momenta, will give different O(e ) terms in eqs. (4.54) and (4.55). However, the 
difference must always be the same; otherwise a 1/e UV pole would be generated from a UV 
convergent integral.) Including the overall factors, we get, 
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Similarly, the iV* double-trace contribution is obtained using the same formula (4.53) with 
P51 replaced by Pn 1 , taken from eqs. (4.47) and (4.49), 
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We notice that, similarly to the numerator factor N^ ad in eq. (4.51), both eqs. (4.56) 
and (4.57) have manifest permutation symmetry. 

Plugging eqs. (4.56) and (4.57) into the full amplitude, including the double-trace part 
of the color factors, the sum over all 24 permutations, and the overall prefactor, we obtain, 
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We have also evaluated the difference between the one-loop outer integrals of hi,h2 and 
J72 and their corresponding MS subtraction terms, and Ij2 h , prior to reduction to 

bubble integrals. That is, we directly integrated eq. (4.50) for general external momenta 
and separately the outer integral of J51 and its subtraction for two special kinematic points. 
In this way we verified the required cancellation of the nonlocal divergent terms ln(— Sy)/e. 

Finally, we add the contribution (4.41) from the graphs 1-50 that have no subdiver- 
gences, in order to obtain the total four-loop divergence (after subtraction of three-loop 
subdivergences): 
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Of course, the double-trace part of the four- loop counterterm must be chosen to cancel 
these poles, 
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corresponding to a nonvanishing divergent coefficient for a counterterm of the schematic 
form, Tr(T> A ~ k F 2 ) Tr(V k F 2 ). (The covariant derivatives may be distributed among the two 
traces in various ways.) It is possible to perform finite shifts of the coefficient of the single- 
trace operator Tr(V 2 F A ) which has a divergent coefficient at three loops. One can trace 
through the effect of such a shift by shifting P^- and P^l' by O(e ) constants. In principle, 



one could remove the 1/e terms in eq. (4.59) that are proportional to (s 2 + t 2 



- it- 



Howe ver, 



one cannot remove the 1/e 2 pole. Nor can one remove the term proportional to (N 2 C3+25Cs), 
because the dependence on color and kinematics is different from the one induced by the 
three-loop counterterm. 

In conclusion, the double-trace terms in the four-point M = 4 sYM amplitude do diverge 
at four loops, saturating the double-trace finiteness bound of D c = 4 + 8/L. 



V. CONCLUSIONS 



In this paper, we recomputed the four-loop four-point amplitudes of M = 4 sYM theory 
and M = 8 supergravity, first obtained in refs. [10, 11]. By exploiting the conjectured duality 
between color and kinematics [29, 30] we found greatly simplified representations. It also 
allowed us to find the form of the complete amplitude, including nonplanar contributions, 
using only planar cut information as input. We confirmed the correctness of the construction 
by comparing the unitarity cuts of the new expressions to the cuts of the earlier forms [10, 11]. 
This provides new nontrivial evidence in favor of the duality conjecture and the associated 
gravity double-copy property. An important advantage of the current construction is that 
once the sYM amplitude has been arranged into a duality-satisfying form, the construction 
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of the corresponding supergravity integrand is trivial: one simply replaces the color factors 
with kinematic numerator factors in each graph. 

The new form of the four-loop four-point amplitude of M = 8 has an important advantage 
over the previous one [10], because no integral displays a worse UV power counting than 
the complete amplitude. This feature greatly simplifies the extraction of the UV divergence 
of the four-loop M = 8 supergravity amplitude in the critical dimension D c = 11/2, corre- 
sponding to the lowest dimension where both M = 4 sYM theory and M = 8 supergravity 
first diverge at four loops. To carry out the required integration we used techniques similar 
to those described in refs. [9-11]. Our results prove that the four-loop four-point amplitude 
of M = 8 supergravity does indeed diverge in the same critical dimension as the correspond- 
ing amplitude of M = 4 super- Yang- Mills theory. Thus, the M = 8 supergravity finiteness 
bound [10] is, in fact, saturated at four loops. The amplitude divergence in D c = 11/2 
means that the T> 8 R A supergravity counterterm has a nonzero divergent coefficient, in much 
the same way as the Tr(D 2 F 4 ) counterterm of M = 4 sYM has a nonvanishing divergent co- 
efficient in this dimension. Moreover, we found that the four-loop finiteness bound [11, 64], 
D < 6, for the double-color-trace terms of M = 4 sYM theory is also saturated. In other 
words, the corresponding D = 6 double-trace counterterms Ti(V i ~ k F 2 ) Ti(V k F 2 ) are also 
present with non-vanishing coefficients. 

More generally, the duality between color and kinematics offers the promise of carry- 
ing advances from the planar sector of gauge theory to the nonplanar sector and then to 
gravity theories. Its underlying origin is, however, still poorly understood; recent progress 
suggests that, at least in the self-dual case [46], underlying it is an infinite-dimensional Lie 
algebra of area-preserving diffeomorphisms. Progress has also been made in finding explicit 
representations of tree amplitudes that manifestly satisfy the duality [47]. It would be inter- 
esting and very useful to devise effective rules that would generate directly duality-satisfying 
representations for loop amplitudes, thus eliminating the need to solve the system of du- 
ality constraints on a case-by-case basis. A step towards finding a Lagrangian with the 
desired properties has been given in ref . [45] . It would also be interesting to explore whether 
the color-kinematic duality extends beyond weak-coupling perturbation theory as well as 
whether the existence of such a duality has practical consequences after carrying out the 
loop-momentum integrations. 

Explicit calculations often lead to surprises. The results described here are no different. 
In particular, in the critical dimension D — 11/2 we found that, after reducing the integrals 
containing UV divergences to a basis of vacuum integrals encoding the numerical factors in 
front of the divergent operator, the UV divergence is given by exactly the same combination 
of basis integrals as found in the single-trace 1/A^-suppressed terms of M = 4 sYM theory. 
It seems unlikely that this is accidental because similar behavior is found at lower loops. 
It would obviously be important to understand the origin of this curious connection and 
implications it may have at higher loops on UV divergences. Another interesting property 
is the existence of strong cancellations between the contributions of various graphs to the 
UV divergence in the critical dimension. This suggests that different integral contributions 
may be related to each other by a hidden symmetry. 

In summary, the duality between color and kinematics offers a powerful means for stream- 
lining the construction of multiloop amplitudes, carrying advances in the planar sector to 
the nonplanar sector. It allowed us to express the numerators of the four-loop four-point 
amplitudes of M = 4 sYM theory and M = 8 supergravity in terms of the numerators of two 
planar graphs. Using this simplified form, in the critical dimension D c = 11/2, we found 
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a surprising coincidence between the UV divergences of M = 8 supergravity and those of 
subleading color single-trace terms of M = 4 sYM theory. This hints at further new relations 
between gauge and gravity theories to be unraveled and that further surprises await us at 
five and higher loops. We look forward to using the tools described in this paper to further 
explore the multiloop structure of gauge and gravity amplitudes and to unravel their UV 
properties. 
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Appendix A: Useful numerator functional equations 

In this appendix we list a set of numerator equations that determine the four-loop four- 
point M = 4 sYM amplitude (up to snail contributions), starting from the two planar master 
graphs 18 and 28 in fig. 4. These equations follow directly from the dual Jacobi relations. 
However, to make the equations more convenient for generating numerator factors from the 
our two planar master numerators, we performed various simplifications which follow from 
the M = 4 sYM auxiliary constraints described in section II B. In particular, we use the 
two-term relations (see eq. (A5)) which rely on the no one-loop triangle subgraph constraint 
to eliminate numerators appearing in other dual Jacobi relations. We also simplified the 
functional arguments of the numerators using the auxiliary constraint that numerators are 
independent of the loop momenta of one-loop box subgraphs. For example, instead of the 
dual relation, 

N5o(ki,k2,k3,ls,lQ, 1-1,1$) = N 2 %(k 2 , ki, k A , 1$, k% — I7, k 2 — 1%, h) 

- N 28 (k 1 , k 2 , k 3 , l 6 , k 4 - l 8 , k x - l 5 , l 7 ) , (Al) 

we simplify this to 

^5o(^i> k 2 , ^5) ^6) hi h) — N 2 s(k2, k±, k^, Z5, k 3 — I7J7, Is) 

- N 28 (ki, k 2 , k 3 , l 6 , k A - l 8 , h, h) , (A2) 
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using the fact that A^s is, in fact, independent of the values of the last two arguments since 
these momenta are those of one-loop box subgraphs (see fig. 4). In this sense, the last two 
arguments of N 2 s are effectively placeholders, and can be assigned any value without altering 
the numerators. These simplifications, however, imply that the equations given below are 
specific to M = 4 sYM theory and will not hold for corresponding numerators of amplitudes 
of theories with fewer supersymmetries. They are also not in direct correspondence with 
the color-Jacobi equations, because numerators of graphs with triangle subgraphs are set to 
zero, although corresponding color factors are nonvanishing. 

On the left-hand side of each duality equation, for simplicity, we will suppress the canon- 
ical arguments, which are the three external momenta and the four independent loop mo- 
menta following the graph labels in figs. 5-11, i.e. 



and we take k± = —k± — k 2 — k% throughout. We have ordered the equations so that the 
substitutions that are required to express the given numerators in terms of the two master 
numerators always come from previous equations in the list. With the above notation, the 




(A3) 
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required equations are 

N 58 = N ls (h, k 2 , h, k 2 - k, k, k, k) ~ N 18 (k 2 , h, k 3 , k x - l 6 , l 5 , l 7 , l 8 ) , 

N 33 = N 28 (h, k 3 , k 2 , k 3 - l 5 , k 2 -l 6 + l 7 , l 7 , l 8 ) - N 18 (ki, k 2 , k 3 , k 2 - Iq, k 3 - k, l 7 , l 8 ) , 

N50 = N 28 {k 2 , ki, k 4 , Z5, k 3 — l 7 , l 7 , l 8 ) — N 28 (ki, k 2 , k 3 , Iq, k 4 — l 8 , l 7 , l 8 ) , 

Nq = ~~ N 33 (kx, k 2 , ki, l 7 , l§ — Iq, k\ — Iq, /§) — N 33 (k 2 , k±, k 4 , l 7 , Iq, k 2 — l 5 + Iq, l 8 ) , 
N 14 = -N 33 (k 3 , k 2 , h, k, -l 5 - l 7 , k 3 -l 7 + l 8 , l 6 ) - N 33 (k 3 , k 2 , h, l 5 , k 2 + l 7 , l 7 - l 8 , l 8 ) , 
N 24 = —N 28 (ki, k 2 , k 3 , 1$ — l 7 , —Iq, l 7 , l 8 ) — N 33 (ki, k 2 , k 4 , —Iq, — 1 7 , —l 5 , l 8 ) , 
iV32 = -N 28 (k 4 , h, ki, l 7 , k 3 - l 5 , l 7 , l 8 ) - N 33 (k 2 , k x , k 3 , l 5 , l 6 , k 2 + l 5 + l 6 - l 7 , l 8 ) , 
^48 = N 28 (k 3 , h, h, l 8 , k 2 - l 5 , l 7 , Z 8 ) - N m (kx, h, h, k 3 - l 6 , k 2 - l 5 , l 7 , l 8 ) , 
jV 49 = _]V 3 3(A;i, k 2 , k 3 , k 3 - l 8 , k 2 - l 5 , -l 7 , l 8 ) - N 33 (k 4 , h, k 2 , l 5 , -l 7 , Iq, l 8 ) , 
^66 = N 58 \ki, k 2 , k 4 , h — k 3 — l 6 , l 6 , l 7 , l 8 ) — N 58 (k x , k 2 , k 3 , k 3 + Iq, Iq, l 7 , l 8 ) , 

Ni = -N e (ki, k 2 , k 3 , Iq, l 5 , l 7 , l 8 ) - N 6 (ki, k 2 , k A , Iq, Z 5 , l 7 , l 8 ) , 
Nq 8 = A r 14 (/c 1 , k 2 , k 3 , k\ — Z5, — Iq, —l 7 , — /§) — N 14 (k x , k 2 , k 4 , l 5 — k 2 , —l 7 , —Iq, l 8 ) , 
N 21 = -N u (k 2 , h, k 3 , k, Iq, l 7 , l 8 ) - N 18 (k 2 , k x , k 3 , -l 5 , k x + k 3 + l 5 - Iq, l 7 , l 8 ) , 
N 2 q = N 24 (k 2 , k\, k 3 , —l 5 , —k 4 — Iq — l 7 , l 8 , Iq) — N 24 (k 2 , k\, k 4 , — 1$, l 7 — k 3 , Iq — ki — l 5 — l 8 , Iq) , 
N 27 = -N 18 (k 2 , ki, k 4 , -h, l 7 , l 7 , l 8 ) - N 2 4(k 1} k 2 , k A , / 5 , -k 3 - l 7 - l 8 , k 3 -l 6 + l 7 + l 8 , l 8 ) , 
N 37 = -N 28 (k 2 , fci, k 3 , fci - h, k A + k, h, k) - N A9 (k 2 , k x , k 3 , k x - l 5 , -l 8 , l 7 - k 2 , l 6 ) , 
N 3 9 = N 28 {k 2 , k\, k 3 , —I5 — l 7 , k A + Iq + l 8 , l§, Iq) — N A8 (ki, k 2 , k 3 , l 7 , l 8 , —1$ — l 7 , —Iq — l 8 ) , 
N45 = N A9 (ki, k 2 , k 3 , l 5 — Iq — l 7 — l 8 , k A — Iq, 1$, l 7 ) 

+ N A9 (ki, k 2 , k 4 , k 2 + l 6 + l 7 + l 8 , l 7 , l 5 , k A - Iq) , 
N 38 = N 49 (k 2 , k x , k A , Iq, k 3 + l 5 + l 7 , -/ 5 + Iq, k 4 - l 8 ) 

— N A g(ki, k 2 , k A , 1$ — Iq, k 3 + l 5 + l 7 , —Iq, l 7 + l 8 ) , 

^53 = N§ 8 (ki, k 2 , k 3 , k 3 — l 8 , Iq, l 7 , l 8 ) + Nqq^i, k 2 , fc 4 , l 8 , —k A — l§, l 7 , l 8 ) , 

N 12 = N 18 (k 4 , k 3 , k 2 , Iq, k 2 + l 8 , l 5 , l 7 ) + A^eOs, h, h, -Iq, l 8 , -/ 5 , k) , 

N51 = N 18 (k 3 , k 2 , k x , k ± + k 2 - l 5 , -Iq, l 7 , l 8 ) - N 21 (k 2 , k 3 , k x , k - h - k 2 , -Iq, l 7 , l 8 ) , 

N 63 = N 21 (k u k 2 , k 3 , k 2 - l 5 , ki + k 2 -l 5 - Iq, l 7 , l 8 ) 

— N 21 (k 2 , kx, k 3 , k x - l 5 , h + k 2 -l 5 - Iq, l 7 , l 8 ) , 
N 7Q = N 45 (ki, k 2 , k 3 , k 2 — 1$, k A — l 7 , Iq, —Iq — l 8 ) 

— N^[k\, k 2 , k 3 , 15 — k\, l 7 , k 3 — Iq, k 4 + Z5 — l 7 — l 8 ) , 

N 80 = N 53 (k u k 2 , k 3 , k 3 - l 7 , Iq, l 7 , l 8 ) + N b3 {k x , k 2 , k 3 , l 7 - k 4 , h, k, k) , 
^55 = N 51 (h, k 2 , k 3 , ki + l 5 , Iq, l 7 , l 8 ) - N 51 (h, k 3 , k 2 , k x + h, Iq, l 7 , l 8 ) , 
N 83 = -N 55 (k 3 , k u k 2 , h + k 2 - l 5 , /§, Iq, l 7 ) - N 55 (k 3 , k u k 2 , l 5 - k 3 , Iq, l 7 , l 8 ) . (A4) 

There are also a set of simpler two-term relations whenever one of the three numerators 
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vanishes due to the appearance of a forbidden one-loop triangle subgraph, 
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Plain-text, computer-readable versions of both the original duality relations (using only the 
no one-loop triangle subgraph property) and the simplified ones presented above may be 
found online [38]. Many other functional equations can be obtained from the dual Jacobi 
relations, which we will not list here. Although important, as they provide additional in- 
dependent constraints to the full system, they are not needed to specify the solution once 
the system and master graph numerators have been solved. However, we have confirmed 
that the numerators presented in appendix B automatically satisfy all these remaining dual 
Jacobi relations. 



Appendix B: Explicit numerators for graphs 

In this appendix we give the explicit values A, of the distinct graph numerators in 
the M = 4 sYM amplitude. (The remaining ones are given directly in terms of these via 
eq. (A5).) These values are obtained by taking the numerators of the master graphs (3.14) 
and substituting their values into eq. (A4). We have performed some algebraic simplifica- 
tions to obtain the results collected here. The M = 8 supergravity numerators are squares 
of the Aj, as in eqs. (3.2) and (4.16). 
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The explicit values of the distinct numerators are 
N e = |s 2 (r 45 - r 35 - s) , 

N 12 = |s(s(7l6 - T"26 - T"35 + ^45 + 2t 56 + 2£) 

- 2(4ti 6 t 25 + 4r 15 r 26 + n 5 (r 36 - 3r 46 ) + T35(r 46 - 3r 36 ))) , 
iVu = \{s{9t? s + 9r 2 2 5 + 4i 2 + 8tr 35 + 2r 2 5 + 2r 4 2 5 ) + 8r 25 (u 2 - s 2 ) 

- 5t(r 15 T 35 + 725T45) + n(4t(2ri 5 + Z 2 ) - 5ri 5 r 45 - 5^5X35)) , 
N 18 = ±(6u 2 r 25 + u(2s(5r 25 + 2r 26 ) - r 15 (7r 16 + 6t)) 

+ t(T 1S T 26 - r 25 (r 16 + 7r 26 )) + s(Ar 15 (t - r 26 ) + 6r 36 (r 3 5 - r 45 ) 

- ri 6 (4t + 5r 25 ) - r 46 (5r 3 5 + r 45 )) + 2s 2 (t + r 26 - r 35 + r 36 + r 56 )) , 

iV 21 = i(t(12r 2 5 - 7r 15 r 16 + r 25 (r 16 - 10r 35 )) + w(t 25 (12t 25 - 8t - 7r 26 ) + r 15 (8t + r 26 - 10r 35 )) 
+ 4/ 2 (m 2 - st) - s(2r 15 (6r 16 + 5r 25 ) + 4w(r 16 - t + r 26 + 2r 45 ) + r 35 (r 36 - 12r 35 - 10r 46 ) 

- r 45 (llr 36 + 6r 46 ) + 12r 25 r 26 ) + 2s 2 (u - r 16 - r 35 + r 36 - r 56 )) , 

^24 = i( s ( 2r 36( 3r 35 + T45) + t- 46 (t 35 + 13t 45 ) - 4u(2r 15 + n 6 - 2r 25 + r 26 ) - llr 15 r 26 ) 

- 2s 2 (2ti 5 - 2r 25 + r 26 - r 36 - r 37 + r 47 - r 56 ) + u(11t w t 25 + ti 5 (7ti 6 + r 26 )) 
+ tr 2 5(12r 16 + 7r 26 )), 

^26 = i( UTl 5( 7r i6 + 12r 26 ) + t(llri 5 r 26 + t 25 (ti 6 + 7r 26 )) + s(16ti 5 Ti 7 

- 4u(2ti 5 + 7i6 - 2r 25 + r 26 ) + r 25 (16r 27 - 11t 16 ) + r 35 (6r 36 - 4r 37 + r 46 - 20r 47 ) 

+ t- 4 5(2t 36 - 20r 37 + 13r 46 - 4r 47 )) + 2s 2 (Yi 7 - 2ri 5 + 2r 25 - r 26 - r 27 + r 36 + r 56 + 2r 57 )) , 
N27 = \(ut 18 (t 25 - 7ri 5 ) + r 28 (t(ri 5 - llr 25 ) - 4wr 25 ) + s(4ti 5 (3ti 7 + r 18 - r 27 ) 

+ 4u(2ti 5 - 2r 25 + r 28 ) + r 45 (5ri 8 + 5r 28 - 16r 37 - 5r 38 ) + r 35 (2r 38 - 16r 47 - 9r 48 ) - 4ri 8 t 

- 4t 25 (ti 7 - 3r 27 )) + 2s 2 (2ri 5 + 2r i7 - 2r 25 + r 28 - r 36 + 2r 37 + r 38 + r 46 + 2r 57 + r 58 )) , 
N 28 = l(s(2r 15 t + r 16 (2t - 5r 25 + r 35 ) + 5r 35 (r 26 + r 36 ) + 2t(2r 46 - r 56 ) - 10wt 25 ) 

- 4s 2 r 25 - 6w(r 46 (t - r 25 + r 45 ) + r 25 r 26 ) - t(r 15 (4r 36 + 5r 46 ) + 5r 25 r 36 )) , 
N32 = l(t{r 25 (r 16 - llr 25 - 12r 26 ) + r 25 r 35 - 6r 3 2 5 - ri 5 (r 26 - 4r 45 )) - u(5t 25 t 2G 

+ 7"i5(7r 16 - 5r 25 + 5r 35 )) + s(t 15 (5t 16 - At) + r 16 (8t + r 25 ) + r 35 r 36 + 5r 45 (r 25 - r 35 - 2r 36 ) 

- r 46 (llr 35 + 6r 45 ) + 2u(6r 25 + 4r 26 - 4r 35 - 2/ 2 )) + 2s 2 (3r 25 - r 35 - r 36 - 3r 46 + r 56 )) , 
N33 = j(s 2 (At 17 - 2(4r 26 + r 35 + 2r 36 )) - 6n 2 r 35 + u((4r 16 + 5r 26 )r 45 

- ri 7 (llr 25 + 7r 35 + 6r 45 ) + r 35 (llr 37 - 5r 46 )) - t(5ri 7 r 25 + 6r 15 r 2e + (6r 26 - 5r i7 - 4r 27 )r 45 
+ 7- 35 (7r 26 - llr 36 + 5r 47 )) + s(ti 5 (5t 26 + 4r 46 ) - 5r 35 (ri 6 + r 27 ) + 4r 25 r 47 + 2w(5r 17 + 2r 25 

- 5(r 26 + r 35 ) + r 56 ) + 2t(r 16 - n 5 - r 27 + r 57 ))) , 

N37 = \(u 2 (^r 15 - 2r 27 ) - 2s 2 (2ri 5 + 3r 27 + 4r 36 ) + t(6r 26 r 27 + 5r 27 r 35 - 6r 27 r 36 

+ 6r 36 r 37 - r 25 (5r 36 + 4r 46 ) - t 15 (6t 26 + 5r 46 ) + 4r 35 r 47 ) + n(5r 27 r 46 - 6r 35 r 36 + (5r 27 
+ 4r 37 )r 45 + r 17 (5r 36 + 4r 46 - 2t) + ti 5 (12t 36 + 5r 37 + 6r 46 - 5r 47 ) + 2t(r 57 - r 25 )) 
+ s(6r 25 r 27 - (4r 26 + 5r 36 )r 45 - 5(r 26 + r 46 )r 47 - T 16 (4t - 5r 27 + r 47 ) - t 15 (11t 26 + 12r 27 
+ 6r 37 + llr 47 ) + 2t(r 35 - r 56 ) + 2u(t + 2r 26 - 5r 27 - 8r 36 + 2r 37 + r 67 ))) , 
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^38 = l(t(T w T 25 + 6r 2 2 5 - llT25T 26 - T 25 T 35 + llrfg + 7X35X37 + T 37 r45 + r 15 (r 45 - 4r 16 ) 

- - s(5rf 5 - 6r 17 r 25 + 20r 17 r 26 + 6r 25 r27 + 4r 26 r27 + 4ri 6 (r 17 + 5r 27 ) - 9t2 5 t 36 

- Ht"35T"36 - 16r 36 r 37 - 5r 37 r 4 5 + 10r 35 r4 6 + 5t 45 t 46 - 4(r 35 + 4r 46 )r 47 + ri 5 (r 37 + r 47 

- 4r 27 + 10r 35 - 9r 36 )) + u(t 45 (7t 47 - 5(2r 35 + r 45 )) - 4r25r 26 + r 15 (r 26 - llr 16 )) 
+ 2s(4tr 37 + u(4(n 6 + r 25 - r 26 - r 35 + r 47 ) + 2lf) + s(2r 16 - 2r 17 + 2r 25 - 2r 26 

- 2r 27 - 2r 35 + r 56 - r 57 + 2r 67 ))) , 

^39 = z(t(r 15 (12r w + 6r 26 + r 36 ) - 10r 38 r 45 + t 27 (5t2 6 + 12(r 36 + r 46 )) + r 35 (r 48 - 7r 38 )) 

- w(t 17 (7ti 6 + r 26 ) + llr 16 r 27 - r 2 5(4r 36 + 5r 46 ) + r 45 (9r 38 + 6r 46 + 7r 48 )) + s(4r 17 r 18 
+ 11ti7T 2 6 + 20ti 8 t 27 + 20ri 7 r 28 + 4t 27 T2 8 + t 15 (5t w - llr 18 + 10r 26 + 5r 28 ) - 4ri 8 r 35 
+ 4r 26 r 3 5 - 6T35r 38 - 16r 37 r 38 + 2t(r 25 + r 36 + 4r 38 ) + 4r 28 r 45 - 5r 25 r 46 - r 37 r 46 

- Il7- 45 r 46 - 13r 46 r 47 - 2r 36 (3r 37 + r 47 ) + 2(3r 35 - r 45 - 8r 47 )r 48 
+ 2m(3ti 5 + 4ri 7 - 4r 27 + 3r 46 + 4r 48 + r 56 )) 

- 2s 2 {u - Ti 5 - 2ri 7 + 2r 18 + r 26 + 2r 27 + 2r 28 + r 35 - r 46 + r 58 + r 67 + 2r 78 )) , 
N 45 = |(8w 2 r 15 + t(ri 5 (ri 6 + 7r 17 + r 26 + r 27 - 5r 35 + 5r 45 ) - (r 16 + r 17 + 8t)r 25 ) 

+ s(r 15 (12r 17 - 18r 25 + llr 26 + r 27 ) + 2(5r 16 r 25 - 2r 17 t - 2t(6r 25 + r 26 ) + 3r 2 5(r26 + 2r 27 )) 
+ llr 3 2 5 + r 45 (6(r 45 - 2r 46 - r 47 ) - r 36 - llr 37 ) + r 35 (r 37 - 5r 36 + 13r 45 - 10r 47 )) 
+ m(4s(ti 6 - i + r 27 + 2r 35 ) - 6r 15 r 16 + r 2 5(7r 26 + 7r 27 - 5r 35 + 5r 45 ) + 4fi 2 ) 
+ 2s 2 (r 16 + r 27 + 5r 35 + r 36 + r 37 + r 38 + 3r 45 - r 48 - r 56 + r 57 )) , 
iV 48 = f (s(ti 8 (5t 25 + 4r 35 ) + r 36 (r 15 - 5r 27 - 4r 45 ) - (t 17 - 5r 25 )r 46 + 5r 35 r 48 ) - 2sut 

- 5ri 5 ri 8 t - t(5r 17 r 26 + 7r 35 r 36 + t 25 (11t 26 - 5r 38 ) + r 15 (5r 28 + r 38 ) + 4r 37 r 46 

+ t 36 t 47 - 4ti 6 (t 45 + r 47 )) + u(4ti 5 T26 + 6r 25 r2 8 + t 16 (7t 17 + r 25 - 4r 35 ) + r 36 (llr 37 
+ r 45 ) - At 2& (t 17 + r 47 ) + 6(r 45 - r 25 )r 48 ) + 2(s 2 (r 25 - r 35 + 2r 36 ) + 3ut(r 17 + r 48 ) 
+ s(u(3t 25 + r 26 + 4r 36 - r 56 ) + t(r 15 - r le - 2r 17 - r 18 - r 47 - 2r 48 + r 58 + r 67 )))) , 
iV 49 = i(2s 2 t - t(llr 3 5r 38 + 2r 15 (2r 36 + 5r 46 ) + t 25 (11t 28 - 4r 37 + 4r 38 + r 47 ) 

+ 4r 27 r 48 + 5ti 7 (t 3 5 + r 48 )) + s(ti 5 (11ti 7 + 5r 38 ) - 5ri 7 r 28 - 2ti 6 (5t 2 5 + 2r 35 ) + r 45 (5r 46 

- 4(7i8 + T"37)) + 5r 38 (r 27 - r 47 ) - 4r 28 r 47 + r 2 5(7r 27 + 5r 48 - 5r 26 )) + M(r 27 (llr 38 + 4r 45 ) 

- t 25 (5t w - 5r 18 + llr 26 + r 37 ) - t 17 (6t 18 - 12r 38 + r 45 ) + (7r 25 + 6r 35 )r 46 + r 38 (5r 45 + 6r 47 ) 
+ t-i 5 (t 46 - 4(r 47 + r 48 ))) + 2(s 2 (r 25 - 4r 17 - r 35 + 2r 38 ) + u(4ur 17 + t(2r 27 - 2r 15 

- 2t 25 - 3r 46 + r 57 )) + s(ri 6 i + u(t 28 + 4r 38 - r 58 ) - i(ri 8 + 2r 37 - 2r 46 + r 56 + r 78 )))) , 
^50 = |(*(t25(5t 3 7 + 4r 47 ) - r 38 (4ri 6 + 5r 26 ) - 5t 15 t 27 ) + s(r 45 (4r 27 + 5r 37 ) 

- ri 8 (5r 26 + 4r 36 ) - 5r 15 r 47 + 5(ri 6 - r 36 )r 48 ) + M(6r 3 5r 37 - 6r 26 r 28 + r 15 (llr 17 - r 37 ) 
+ (llr 16 + 6(2r 26 + r 36 ))r 48 ) + 2t(3u(r 15 - r 26 + r 37 - r 48 ) + s(r 18 - 2r 15 + 2r 26 - r 27 
+ r 36 - 2r 37 - r 45 + 2r 48 + r 57 - r 68 ))) , 

N51 = |(4t 2 (ri 5 + ^25) - 6r 2 5 t - 4st(s - u) - 6ut% 5 + T 35 (5tr 25 - 6sr 35 ) 

+ 5t 15 (st 25 + UT35) + 14s 2 r 45 - s(6uti 5 + t(n 5 + 6r 25 + 13r 35 - 2/ 2 )) - 2u 2 Z 2 ) , 
j\r 53 = 8s(tr 35 + «r 45 - sr 25 ) , 
^55 = lt(t(r 25 - 8r 15 + 5r 45 ) + w(9r 45 - 17r 15 )) , 
N 58 = s(2w(r 45 - 3r 35 ) - s(u - t + 4r 25 + 5r 35 + r 45 )) , 
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N 63 = -H5tr 35 + M ( 12r 36 + 5r 45 - 4r 46 ) + s(2t - 2u + 2r 25 + 2r 46 + 8r 26 + 10r 36 - 7r 15 )) , 
iV 66 = s(4t(r 35 - 2r 36 ) + 2u(r 35 + 3r 45 - 4r 46 ) - s(6u + t 15 - 6t + 5r 25 - 8r 26 )) , 
AW = |s(s(2(r 15 + 1 - u) - 7r 25 ) + 5ur 35 + 5tr 45 ) , 
N 79 = §s(5*t 35 + w(5r 45 + 4r 46 - 4r 37 + 12r 47 - 12r 36 ) 

+ s(2(n 6 + r 25 - r 27 + 3r 17 - 3r 26 + 4r 47 - 4r 36 ) - 7r 15 )) , 
N 80 = 16s 2 (u - t) , 

N 83 = -%kls(u-t). (Bl) 

The other numerators are given directly in terms of these in eq. (A5). For Ag 3 it should 
be understood that k\ — > only after canceling the propagator. Alternatively, we can 
rewrite the snail contributions in terms of the numerators of the graphs in fig. 12, we have 

AW = AW = AW = ~s(u - t) . (B2) 

Plain-text, computer-readable versions of these expressions may be found online [38]. 

It is interesting to note that AW can be used as a non-planar master graph numerator, 
as discussed in section III B. This implies that the single numerator AW contains the same 
amplitude-specific information as the two planar master numerators AW and A r 28 combined. 

Appendix C: Vacuum integrals and their expression in terms of master integrals 

In the first column of table I we give the numerators of the vacuum integrals in figs. 18-20, 
as they appear in the expression for defined in eqs. (4.20) and (4.21). Each integral can 
be reduced to a linear combination of the three master integrals Vi, Vi and Vs. The second, 
third and fourth columns of the table provide the coefficients of Vi, Vz and V$, respectively, 
after this reduction. To obtain the coefficient of each vacuum integral Vi in the final formula 
for the four-loop UV divergence in eq. (4.37), we simply sum the numbers in each column 
labeled by a Vi, to obtain 

~ 23 

V (4) = — (Vi + 2V 2 + V 8 ) . (CI) 
Inserting this value into eq. (4.20) yields the final result (4.37). 
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TABLE I: Effective numerators for the vacuum integrals If entering the UV pole of the four- 
loop J\f = 8 supergravity amplitude, and the coefficients arising from writing them as a linear 
combination of basis vacuum integrals V\, V2, and V$. 
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